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Abstract 


A  robust  multivariable  controller  design  technique  is 
outlined.  The  problems  encountered  when  applying  this 
technique  to  a  nonlinear  system  are  analyzed.  To  overcome 
these  problems,  it  is  proposed  to  apply  the  robust 
controller  in  a  cascaded  scheme  with  simple  multi  loop 
proportional  or  propor t i ona 1 -p 1 us-der i vat i ve  control  in  the 
slave  loop.  The  performance  of  this  control  arrangement  is 
evaluated  on  a  nonlinear  distillation  column  simulation 
model  and  compared  to  conventional  PID  control. 

To  apply  the  robust  multivariable  controller  to  a  pilot 
plant  distillation  column,  a  control  system  utilizing  an 
LSI- 11  microcomputer  has  been  developed.  The  system  software 
and  hardware,  designed  to  provide  flexibility  and  a 
comprehensive  operator  interface,  are  described  in  detail. 
The  modifications  necessary  to  adapt  the  system  to  different 
control  algorithms  are  indicated. 

The  newly  developed  control  system  is  used  to  apply 
multivariable  robust  control  and  multi  loop  PID  control  to  a 
pilot  plant  distillation  column.  The  robust  controller  is 
used  both  with  proportional  and  proport ional -plus-derivative 
control  in  the  slave  loop. 

Comparison  of  the  results  shows  that  the  robust 
controller  generally  performs  better  than  conventional 
multi  loop  PID  control. 
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1.  Introduction 


There  is  a  considerable  interest  in  the  design  of 
controllers  for  multivariable  systems  [2,3].  Most  of  the 
modern  control  synthesis  techniques  require  a  mathematical 
model  of  the  process  that  is  to  be  controlled.  In  many 
instances  the  modelling  of  complex,  often  poorly  understood 
processes  is  a  substantial  and  crucial  task  and  by  no  means 
routine.  A  promising  alternative  for  the  control  of  such 
processes  is  to  employ  design  methods  that  require  only 
plant  input-output  data  and  not  an  explicit  mathematical 
model  of  the  process.  Self- tuning  methods  that  have  arisen 
largely  out  of  the  work  of  Astrom  and  co-workers  [4-6]  have 
been  used  successfully  to  control  a  number  of  industrial 
processes.  A  more  recent  technique  introduced  by  Davison  [1] 
relies  only  on  plant  input-output  data  and  is  based  on  the 
notion  of  compensator  identification  as  opposed  to  plant 
identification  [7-10].  It  is  to  be  noted  that  Davison's 
controller  and  Francis  and  Wonham' s  [11]  controller  based  on 
the  internal  model  principle  are  very  similar  in  structure 
in  that  both  require  an  internal  model.  A  good  comparison  of 
the  two  methodologies  is  described  by  Kwatny  and  Kalnitsky 
[12]. 

Two  limitations  apply  to  Davison's  technique:  the  plant 
char acter i s t i cs  should  be  such  that  it  can  be  described  by  a 
linear,  t ime- i nvar i ant  model;  and  it  is  also  required  to  be 
open- loop  asymptotically  stable.  However,  unlike  the  self 
tuning  methods,  Davison's  robust  feedback- feedforward 
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technique  does  not  require  any  assumptions  regarding  the 
plant  model  order  ,  minimum  phase  characteristics,  etc. 
Secondly,  the  technique  guarantees  a  robust  (with  respect  to 
small  parameter  variations),  closed  loop,  asymptotically 
stable,  multivariable  system.  In  contrast,  convergence  and 
stability  properties  of  even  single- loop  systems  under 
self-tuning  control  are  difficult  to  analyze  or  guarantee 
[13]. 

The  availability  of  sophisticated  computing  hardware 
and  the  development  of  complex  control  strategies  have 
triggered  a  rapid  development  of  computer  control  systems 
[19,20,22,23].  The  tendency  is  away  from  single  maxi 
computers  to  hierarchial  networks  of  mini  computers  as  these 
provide  more  flexibility  and  reliability  combined  with  more 
than  enough  computing  power.  The  detailed  development  of  one 
"node"  in  a  control  network  is  necessary  for  the 
implementation  of  complex  process  control  algorithms. 

The  purpose  of  this  thesis  is  twofold:  the  evaluation 
of  Davison's  robust  control  technique  on  a  realistic  plant 
and  the  development  of  the  tools  to  do  so. 

This  thesis  is  organized  as  follows.  In  chapter  2  the 
theory  of  robust  multivariable  control  is  outlined.  The 
difficulties  posed  by  nonlinear  systems  are  indicated  and  a 
simple  but  elegant  way  of  solving  those  problems  is 
proposed.  The  chapter  is  concluded  with  an  extensive 
evaluation  of  the  controller  on  a  simulation  program. 

Chapter  3  describes  the  development  of  a  control  system 


; 
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which  involves  putting  together  the  necessary  hardware  and 
developing  the  software.  Although  the  system  was  developed 
with  the  department's  pilot  plant  distillation  column  in 
mind,  it  can  easily  be  adapted  to  other  systems. 

Chapter  4  is  the  synthesis  of  chapters  2  and  3:  the 
application  of  the  robust  control  technique  on  the 
distillation  column.  An  extensive  series  of  experiments 
forms  the  basis  for  the  controller  evaluation. 

The  last  chapter  contains  the  conclusions  and 
recommendations  for  future  work. 


2.  Theory  of  Robust  Control 


2.1  Introduction 

A  majority  of  the  advanced  control  techniques  deal  with 
idealized  linear  systems  and  are  generally  more  complex  than 
conventional  PID  control  and  require  in  depth  Knowledge  of 
the  plant  that  is  to  be  controlled.  In  comparison  to  this, 
the  design  and  application  of  a  multivariable  tuning 
technique  due  to  Davison  is  simple  and  straightforward . 
However,  like  most  other  methods  it  has  been  developed  for 
idealized  linear  systems,  its  performance  for  nonlinear 
plants  is  unknown  and  a  subject  of  interest. 

At  this  point  we  should  point  out  that  there  seems  to 
be  no  single  definition  of  the  term  "robustness"  in  the 
literature.  Generally  robustness  of  a  control  system 
indicates  the  fact  that  the  system  performance  is  maintained 
in  the  face  of  uncertainties  or  differences  between  plant 
model  and  actual  plant.  Usually  when  we  think  of  robustness, 
we  specify  gain  and  phase  margins  for  single  loop  linear 
systems . 

A  possible  measure  of  robustness  for  a  feedback  system 
is  the  magnitude  of  (otherwise  arbitrary)  per turbat ions 
which  may  be  tolerated  without  instability.  By 
"perturbations"  is  meant  the  deviation  of  the  plant  transfer 
function  model  G(s)  from  the  true  plant  [14]. 

Davison  [1]  does  not  quantify  robustness,  but  proves 
that  asymptotic  tracking  and  output  regulation  occurs 
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independent  of  input  disturbances  and  plant  parameter 
variations,  as  long  as  the  plant  is  stable  and  is 
describable  by  a  finite  dimensional,  linear,  t ime- i nvar i ant 
state-space  model . 

In  this  chapter  we  will  outline  the  theory  as 
originally  published  by  Davison  [1].  Since  the  theory  was 
developed  for  linear  systems,  we  may  expect  to  encounter 
some  problems  when  applying  the  Robust  Controller  directly 
to  a  real  life  nonlinear  system.  We  will  look  at  the 
character  of  the  difficulties  and  its  causes. 

Then  we  will  make  a  proposal  to  deal  with  the  nonlinear 
system  to  make  its  behaviour  acceptable  for  the  robust 
controller  by  i ncorporat i ng  slave  control  loops  in  the  total 
control  system.  Several  advantages  to  this  approach  will  be 
out  1 i ned . 

Finally,  we  will  design  and  evaluate  a  Robust 
Controller  for  a  nonlinear  distillation  column  simulation 
program. 


2.2  Theory  of  Robust  Control 

2.2.1  Davison's  Robust  Feedforward  -  Robust  Controller 
Consider  the  following  linear  t ime- i nvar i ant  stable 
p 1  ant : 

x  =  Ax  +  Bu  +  Ew 
e  =  Cx  +  Du  +  Fw  -  y 

y ( s )  =  C ( s I - A  f*  Bu(s)  +  C(sl-A)  fEw(s)  +  Du(s) 


(1) 


. 
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K  *vi  t 

where  xe R  is  the  state,  ufcR  is  the  control  input,  ee R 
is  the  error  which  is  the  difference  between  the  output 

y  =  Cx  +  Du  +  Fw 

and  the  specified  reference  input  y^p  ,  where  y  ««  Rr  and 

d  *  * 

w*R  is  an  input  disturbance. 

It  is  assumed  that  the  elements  of  w  and  y^e^  satisfy 

the  same  differential  equation 

w(p)  +  apW^  ^  +  .  +  a2w  +  a,  w  =  0  (2) 

which  has  characteristic  roots  all  contained  in  the  closed 
right  hand  part  of  the  complex  plane  and  that  the  initial 
conditions  of  (2)  are  independent  of  each  other  for  each 
element  of  w  and  y ^  . 

We  will  only  consider  step  changes  in  set  points  and 
disturbance  inputs.  For  this  type  of  inputs  (2)  becomes 

w  =  0  (3) 


with  characteristic  root  %,  =  0  which  is  contained  in  the 
closed  right  hand  part  of  the  complex  plane. 

From  Davison's  more  general  theorems  we  can  derive 
directly: 

Theorem  1 

If  the  disturbance  w  is  measurable  then  there  exists  a 
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feed  forward  controller  for  the  system  (1)  so  that  e(t)  = 

l 

(yre^  -yl-^O  as  t  -*■  04  for  all  constant  disturbances  w*R 
and  constant  reference  inputs  y^  e  R1"  if  and  only  if 

rank  G(0)  =  rank  C(sI-A),B  +  d|  =  r  (4) 

I  s 

Theorem  2 

There  exists  a  robust  controller  for  the  system  (1) 

such  that  e(t)  =  (yv.cJ)  ~  y)"*0  as  t-*  00  for  all  constant 

, ( 

disturbances  w«R  and  all  constant  reference  inputs 
Yref6^  ar)d  the  system  is  asymptotically  stable  if  and  only 
if 


(5) 


rank  G  ( 0 )  =  rank  C(sI-A)HB  +  D  =  r 
Davison's  T,1 -matrix  is  equivalent  to  G  ( 0 )  used  in  frequency 
domain  techniques. 


T he  con trol ler  i dent  if i cat  ion  procedure 

If  the  conditions  of  theorems  1  and  2  are  satisfied,  we 
can  proceed  with  the  controller  identification.  We  will 
consider  a  controller  for  tracking  of  constant  set  points 
and  regulatory  control  under  constant  disturbances,  i.e.  the 
controller  will  be  designed  for  step  changes  in  set  points 
and  disturbances,  not  for  ramps,  parabolas,  etc. 

We  have  to  perform  two  sets  of  experiments  on  the  plant 
to  find  the  Steady-State  Tracking  Gain  Matrix  and  the 
Steady-State  Disturbance  Gain  Matrix: 

Experiment  I  (to  find  T't  ) 


^  1  U  n  >  '»cr  '  •  s  .  ^ ~  1  ^ 


8 


Step  1 :  Apply  an  input  u ( t )  =  u,  ,  where  u,e  Rm ,  u,/  0;  then 
since  the  system  (1)  is  stable,  there  exists  a  constant 
vector  y'  so  that 

1  im  { y ( t )  -  y‘  }  =  0 

-4-*oo 

Measure  y[  . 

In  other  words,  measure  the  steady-state  response  yj  to 
a  step  input  u,  . 

Step  2:  Apply  an  input  u(t)  =  Tj2  ,  where  u2*R  and  u2  is 

linearly  independent  of  u,  ,u3 . u^;  then  there  exists 

a  constant  vector  y^  so  that 

lim  { y  ( t )  -  y^}  =  0 

{;-+  oo 

Measure  y^ • 

In  other  words,  measure  the  steady-state  response  y[  to 
a  step  input  Ti  which  is  independent  of  all  other  step 
inputs  in  experiment  I. 


Step  m:  Apply  an  input  u(t)  =  u*,  where  u^R  and  u^  is 

linearly  independent  of  u,  ,uz . u^  (  ;  then  there 

exists  a  constant  vector  y^  so  that 
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1 im  { y ( t )  -  y^}  =  0 

t~»O0 

Measure  y^. 

Then , 


i  - ( y \  ,yL  t  •  •  • » yin )  (  li  ,  ,u^, . . .  jUj^) 

Experiment  II  (to  find  D',) 

It  is  now  assumed  that  the  disturbance  inputs  w  can  be 
excited.  We  perform  exactly  the  same  procedure  as  in 
Experiment  I,  but  this  time  we  excite  the  disturbance  inputs 
instead  of  the  control  inputs: 

Step  1  :  Apply  an  i  nput  w ( t )  =  w,  ,  where  w,  e  R°* ,  w/  t  0 ;  then 

-  J 

we  find  a  constant  vector  y(  so  that 
1  im  { y  ( t )  -  yj  }  =  0 

1-+OG 

Measure  yj. 


Step  d:  Apply  an  input  w(t)  =  wj  ,  where  wj  c- R  and  wj  is 
linearly  independent  of  w,  , wL  ,  .  .  .  , wj_;  . 

Then  we  find  a  constant  vector  y  so  that 

1 im  { y ( t )  -  y j}  =  0 

{.-+04 


Measure  y^. 


Then 


d|  =  (yi.yi. 


.  ,yj)  (w,  ,w^ . w jf 


and  the  feedforward  -  Robust  controller  is  given  by: 

u  -  [T,'1+yvef  •  I t|  ] +  d ; w  -  t t ;  \+e.  f  (y-y„j>  )dt' 


where  T,  =  tracking  steady  state  gain  matrix 

D[  =  disturbance  steady  state  gain  matrix 
y^f  -set  points 
y  =  plant  outputs 
w  =  disturbance  inputs 
£  =tuning  constant 
+  denotes  pseudo  inverse: 

i  t  1 1+=  i  t  ;  iV.Ti'r' 

for  a  square  matrix  [ T  |  ]*  =  [ T  1  ) 


1 1 


2.2.2  Feedforward  -  Robust  Controllers  and  Frequency  Domain 
Techniques 

In  the  previous  section  we  looked  at  the 
feedforward  -  robust  controller  in  the  time  domain.  In  this 
section  we  try  to  translate  some  results  of  (state-space) 
robust  controller  design  to  frequency  domain  design. 

The  feedforward  -  robust  controller  for  a  stable  plant 
and  constant  disturbances  was  given  by  (6), 

U  =  [  T I  Ty*.^  -  [T,1  ]+D'w  -  [t;  }+*f  (y  -  )dt'  (6) 

o 

As  we  saw,  T |  and  d|  can  be  determined  by  some  simple 
steady  state  experiments.  Since  we  are  interested  in 
frequency  domain  techniques,  we  consider  the  relationship 
between  t|,  D*(  and  the  transfer  function  matrices. 

Figure  1  shows  the  plant  configuration. 

Column  1  of  T  |  contains  the  steady  state  response  of  the 
output  vector  y  to  a  unit  step  in  u,  ,  column  2  the  response 
of  y  to  a  unit  step  in  u2 ,  etc.  Hence  it  is  obvious  that  Tj 
is  the  same  as  G ( 0 )  and  in  the  same  way  D[ =  D(0). 

General ly , 


i  t  -i 

i  =  (-1  ) 


.t-i 


( t -  1  )  I  dstw 


G  ( s ) 


(7) 


■ 
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w(s) 


iL_ 

D(s) 


Figure  1.  Plant  representat i on  in  s-domain 


where 

t  =  1,  2 . qk,  qk  is  the  multiplicity  of  root 

k  =  1,  2,  . ..,  p  ,  p  is  the  number  of  different 
character i st i c  roots  of  (2) 

The  same  type  of  relationship  exists  for  d'  ,  replace  G(s)  by 
D  (  s )  in  ( 7 )  . 

We  consider  constant  incuts  and  disturbances,  so  equation 
(3),  w  =  0,  applies.  This  has  one  root,  X,  =  0  and  from  (7) 
we  find, 


. 
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T 


i 


G(s)l  .  =  G  ( 0 ) 

•  S  -  /lt  =  o 


(8) 


We  can  rewrite  (6)  as 

u(s)  =  [G(  0 )  ]+  {yvtp  (s)  -  D  (  0  )w(  s  I  -  f  (y  ( s )  -yft{  ( s ) }  ( 9 ) 
See  figure  2  for  the  block  diagram. 

The  closed  loop  transfer  function  is 

y(s)  =  G(s)u(s)  +  D(s)w(s) 

=  G(s)G(0)+ {y^  (s)  -  D(0)w(s)  -  £.(y{s) 

"  yh*j>  (  s  )  ) }  +  D  (  s  )w(  s  )  (10) 

y(s)  =  [I  +  |g(s)G(0)+  f'G(s)G(0)+  ( 1  +f)y^(s) 

+  [I  +  fG(s)G(0)+]"(D(s)  -  G(s)G(0)+D(0)  )w( s )  (11) 
When  we  apply  step  inputs  to  y^  and  w,  the  final  value 
theorem  shows 

lim  y(t)  =  lim  sy(s)  =  y vo  (12) 

S ~*0  t  ' 


A  possible  design  procedure  for  a  constant  set  point  and 
disturbances  controller  is: 

1.  Using  DNA  (Direct  Nyquist  Array)  and  Character i st i c 


' :  ■ 
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Figure  2.  Block  diagram  of  controlled  plant 

Locus  Techniques,  design  a  controller  K(s)  for  the  plant 
Q(s),  such  that  the  resultant  closed  loop  system  G(s)  is 
stable  (see  fig.  3)  [15,16] 

2.  Using  DNA ,  obtain  G(0). 

3.  Implement  the  control  scheme  as  shown  in  figure  2. 

4.  Tune  the  controller  by  varying  epsilon. 

We  will  briefly  look  at  a  feedforward  controller  for 
constant  disturbances  and  tracking  set  points.  In  this  case 
the  set  points  satisfy  the  differential  equation 


=  0  (13) 

which  has  characteristic  root  /l,  =  0  with  multiplicity 


2. 
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u(s) 


y  (s) 


Figure  3.  Modified  Plant  G(s) 

For  the  control ler  we  need 
T|  =  G( 0 ) 

,  dG ( s ) 

T,  = -  =  -  G'  (0) 

ds  s  =  0 


(  14) 


(  15) 


and  the  controller  is  given  by 

u(s)  =  -  G  (  0  )+ D  (  0  )w(  s  )  +  G  (  0  )+  (  s  ) 

-  G  (  0  )■"  G7  (0)G(0)  syM|  (s) 

-  G(  0  )+ (  ^  +  ^)  (y(  s  )  -  yhef  (s))  (16) 

See  figure  4  for  the  block  diagram 

We  have  seen  that  the  robust  controller  design 
technique  can  be  useful  in  frequency  domain  design.  However 
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^  y (s) 


Figure  4.  Plant  with  controller  for  set  point  tracking 

if  the  desired  outputs  or  disturbances  follow  higher  order 
differential  equations,  the  design  becomes  more  cumbersome 
because  we  have  to  evaluate  derivatives  of  G(s)  and  G(^k). 
This  is  not  easily  done  by  graphical  methods  since  we 
usually  work  with  G(<°). 


- 
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2.3  Application  of  Robust  Controller  to  a  Nonlinear  System 
2.3.1  Introduction 

In  the  following  discussion,  it  will  be  assumed  that 
the  system  under  consideration  is  2x2,  i.e.  two  inputs  and 
two  outputs.  For  a  square  matrix  [T(‘]+=  [  T  J  ] 

As  we  outlined  in  the  previous  chapters,  we  use  steady  state 
data  for  the  design  of  a  feedback  controller: 


u 


r 


i  -i 

[T  ]  € 


J 


(y 


(17) 


If  we  change  the  sign  but  not  the  magnitude  of  a  control 
input,  only  the  sign  of  the  (steady-state)  plant  output  will 
change,  not  the  magnitude.  When  we  design  a  robust  feedback 
controller  for  a  linear  plant,  both  positive  and  negative 
control  inputs  will  lead  to  the  same  controller  as  we  will 
see. 

In  a  nonlinear  plant,  the  effects  of  positive  and 
negative  inputs  will  often  cause  outputs  to  be  different  in 
magnitude  as  well  as  sign.  In  that  case,  we  will  find  two 
different  feedback  controllers  for  the  same  operating  point 
of  a  nonlinear  plant,  one  for  positive  inputs  (disturbances) 
and  one  for  negative  inputs. 

The  most  serious  difference  in  the  controllers  would  be 
a  difference  in  sign,  i.e.  one  controller  would  perform 


. 
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positive  feedback  and  the  other  negative  feedback.  In  that 
case  the  controllers  would  have  to  be  rejected  for  stability 
reasons.  We  will  investigate  under  what  circumstances  a  sign 
change  will  occur.  The  effect  of  changes  in  magnitude  of  the 
entries  in  the  feedback  controller  matrix  will  not  be 
considered . 

2.3.2  Distillation  Column,  General 

Let  the  inputs  and  outputs  of  a  distillation  column  be 
denoted  by: 

u,  =  reflux  flowrate  y,.  =  top  composition 

u2  =  steam  flowrate  y2.  =  bottom  composition 

All  variables  are  perturbation  variables 

y1(  =  change  in  bottom  composition  resulting  from  a  change 

in  reflux  flowrate. 

In  the  controller,  the  most  important  role  is  played  by  the 
matrix 


-I 


u,  ut 


(  18) 


write  (18)  as 


'  r  ^  ac  ^ 
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(19) 


where  K' s  are  steady-state  gain  parameters. 

If  u,  and  are  positive,  then  yn  and  yZ(  will  be 
positive  (i.e.  an  increase  in  reflux  flowrate  will  cause  an 
increase  in  both  top  and  bottom  composition),  and  y  and 
y  will  be  negative. 

Denote  the  steady-state  gain  parameters  for  positive  inputs 

with  a  superscript  P: 

p  p  p  p 

K|(  ,  Ku  are  positive  and  K(l>  ,  K21  are  negative. 

If  u,  and  u2  are  negative,  y1(  and  yu  will  be  negative 

and  yn  and  y  wi 1 1  be  positive,  but 

M  w  N  N 

K(j  ,  K"  are  positive  and  K/l  ,  K1Z  are  negative  . 

Ergo,  the  entries  in  the  steady-state  gain  matrix  have  the 
same  signs  for  positive  and  negative  inputs.  But  for  control 
we  use  the  inverse  of  the  steady-state  gain  matrix. 

The  entries  in  the  inverse  of  the  Kp  matrix  will  have 
opposite  signs  of  the  entries  in  the  K*  matrix  if 
sign(  det  KP )  i  sign(det  K*) 


We  will  now  investigate  what  the  physical  meaning  of  a  sign 
change  is. 

p 

Consider  det  K  : 


det  K 


(20) 
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From  steady-state  data,  10%  change  in  inputs  we  find, 


det 

Kp 

>  o, 

(21  ) 

KP 

ii 

K« 

- 

Kj,  >  0  , 

(22) 

< 

Ki 

>  < 

kT, 

(23) 

Taking  the  signs  of  K' s  into  account,  we  can  write  (23)  as 

<  k  rt 

— r  <  (9 


or  ,  r(p  <  r£  (25) 


The  parameter  r  is  a  gain  ratio  and  can  be  considered  as  a 
measure  of  selectivity: 

rp  >  1  means  control  input  1  affects  the  top  product 
composition  stronger  than  the  bottom  product 
compos i t i on , 

rp  =  1  means  control  input  1  affects  top  and  bottom  product 
composition  equally  strong, 

rj*  <  1  means  control  input  1  affects  the  bottom  product 
composition  stronger  than  the  top  product  composition. 
Expression  (25)  in  words  is: 

for  positive  control  inputs ,  the  select iv ity  of  the  steam 
flow rate  ( input  2)  for  the  top  product  composition  is 
stronger  than  the  select iv ity  of  the  reflux  flowrate  for  the 
top  product  composition. 


" 
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For  negative  inputs  we  find 


det  K 


N 


(26) 


hence 


(27) 


N 


>  r 


N 


or  ,  r 


(28) 


Expression  (28)  in  words  is: 

for  negative  control  inputs ,  the  select iv ity  of  the  reflux 
flowrate  for  the  top  product  compos  it  ion  is  stronger  than 
the  selectivity  of  the  steam  flowrate  for  the  top  product 
compos i t ion. 

We  see  that  the  emphasis  in  the  response  has  changed, 
now  the  most  effective  way  of  influencing  the  top  product 
composition  with  as  little  effect  as  possible  on  the  bottom 
product  composition  appears  to  be  a  change  in  reflux 
f lowr ate . 

Somewhere  between  the  situations  det  K  negative  and  det 
K  positive  we  expect  the  case  of  det  K  =  0.  For  a  linear 
system  this  would  mean  uncontrollability  since  the  outputs 
would  no  longer  be  independent,  i.e.  y,  would  be  a  multiple 
of  y2  for  every  combination  of  inputs.  Flowever ,  in  our 
nonlinear  situation,  det  K  is  not  independent  of  the  sign 
and  magnitude  of  the  inputs  and  we  cannot  speak  of 
"dependent  outputs  for  every  combination  of  inputs". 
Therefore,  in  this  case  the  occurence  of  det  K  =  0  has  only 
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mathematical  significance,  indicating  a  point  where 
Davison's  synthesis  technique  is  not  applicable. 

Simulation  and  experimental  data  elsewhere  in  this 
thesis  will  show  that  we  indeed  get  a  sign  change  in 
controllers  designed  with  Davison's  method  if  we  use 
negative  instead  of  positive  inputs  and  disturbances  in  the 
plant  experiments. 


2.3.3  Distillation  Column,  Numerical  Example 

Some  open  loop  experiments  were  performed  on  a 
nonlinear  distillation  column  simulation  program  (see 
section  2.5.2)  to  illustrate  the  effect  of  non  1 i near i t i es . 
The  inputs  in  this  case  were  reflux  flow  rate  and  steam  flow 
rate,  the  outputs  were  top  product  composition  and  bottom 
product  composition. 

It  was  found  that  for  positive  inputs: 


0.788 

-2.765 

’l  .21 

0 

-a 

4.319 

-3.981. 

0 

1  .458 

) 

0.651 

-1 .896" 

1 

=  I  Kp  ] 

-i 

3.569 

-2.731 

det  KP  =  4.99 

rf  =  0.182 

r,p  =  0.823 
2 
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-0.547  0 . 380" 

.-0.715  0 . 1 3 1  _ 


and  for  negative  inputs: 


jf- 1.291  0.991' 

\L-2 . 335  12.456 


-1.21  0  " 
.  0  -1.458 


-  1 


1 .067  -0.680' 
1 .930  -8.543 


-1 


det  K*  =  -7.80 


r,v  =  0.553 
r  £  =  0.080 


[t|  P 


1.095  -0.087 

0.247  -0.137 


As  we  can  see,  the  signs  in  the  [ T *  ]  1  matrices  are 
opposite,  which  means  that  at  least  one  of  the  matrices 
would  produce  an  unstable  controller  (positive  feedback). 
Also  the  magnitude  of  the  entries  is  very  different  in  both 
matrices.  We  may  conclude  that  the  robust  controller  design 
technique  simply  does  not  work  with  this  plant. 
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2.4  Modifications  for  a  Nonlinear  System 

As  we  saw  in  previous  sections,  the  nonlinear  behaviour 
of  plants  like  the  distillation  column  makes  direct 
application  of  the  robust  controller  impossible  because  the 
controller  design  gives  ambiguous  results. 

We  suggest  a  rather  simple  but  practically  useful 
modification  to  extend  the  applicability  of  the  robust 
design  technique  to  a  larger  class  of  processes. 

Utilization  of  simple  multi  loop  proportional  or 
proportional  +  derivative  control  for  a  typical  plant 
achieves  two  very  important  plant  character i s t i cs .  With  the 
correct  choice  of  gains  the  plant  can  be  made  asymptotically 
stable  and  secondly  the  assumption  of  linearity  is  valid 
over  a  wider  range  of  operating  conditions  than  the  original 
plant,  due  to  the  effects  of  K  in  the  forward  path  and  the 
unity  gain  feedback  loop.  The  robust  controller  can  be 
synthesized  in  such  a  cascaded  scheme  and  implemented  as  a 
master  loop  as  shown  in  figure  5.  The  slave  or  inner  loop 
consists  of  the  original  plant  with  multi  loop  proportional 
feedback  gain  K,  or  proportional  plus  derivative  feedback. 

It  is  not  possible  to  have  multi  loop  PI  control  in  the  inner 
or  slave  loop  because  this  would  eliminate  steady-state 
offset  in  the  outputs  and  result  in  [  T*  ]  =  0.,  suggesting 
non-existence  of  the  robust  control ler  [1]. 

The  above  suggested  cascaded  scheme  has  three 
advantages.  Firstly,  the  proportional  feedback  action  may 
make  the  plant  appear  more  linear,  i.e.,  the  excursions  due 


■ 
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Figure  5.  P-Robust  controlled  plant 

to  set  point  interactions  and  disturbances  are  reduced.  The 
linearity  assumption  can  now  be  considered  to  be  valid  and 
Theorems  I  and  II  of  Davison  and  the  resulting  algorithms 
can  be  applied  to  the  augmented  'linear'  plant.  Secondly,  in 
the  suggested  configuration  it  is  possible  to  design  the 
robust  controller  for  an  industrial  unit  without  having  to 
do  experiments  I  and  II  on  the  open- loop  plant.  The  normal 
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steady-state  input-output  data  would  be  sufficient  for  the 
controller  design.  Thirdly,  the  plant  inherently  has  analog 
backup  in  the  form  of  proportional  (or  PD)  feedback  control 
in  the  event  of  a  computer  failure.  The  two  latter  features 
are  very  important  for  field  implementation  of  advanced 
control  strategies  on  industrial  plant  units. 

For  a  plant  described  by  eqn .  (1),  the  augmented  plant 
with  controller,  u  =  Ku  -  Ky,  becomes 

x  =  (A  -  BKC ) x  +  BKu  +  Ew 

y  =  Cx  (29) 

u  is  the  set  point  for  the  slave  loops 


Thus  the  modified  plant  is  char acter i zed  by  matrices  (A-BKC) 
=  A,  BK  =  B,  and  E  =  E. 

In  case  of  PD  slave  control,  the  augmented  plant  will 
be  slightly  different.  The  control  law  is: 


u  =  -  Kpy  -  K0y  +  Kp  u 

=  -  KpCx  -  K0Cx  +  Kpu  (30) 

where  Kp  is  the  proportional  gain  in  the  slave  loop  and 
K0  is  the  derivative  gain. 

Then , 


x 

(I 


Ax  -  BKpCx  -  BKdCx  +  BKp  u  +  Ew 
+  BK0  C ) x  =  (A  -  BKpC ) x  +  BKp  u  +  Ew 
[I  +  BK0  C  ]“'  (  A  -  BKp  C  )  x  +  [I  +  BKt>C]~,BKpu 


(31  ) 


x 


iq  *  : 

' .  • 
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+  [I  +  BK^Cf'Ew  (32) 

Rewrite  equation  (32)  as: 

x  =  Ax  +  Bu  +  Ew  (33) 

To  guarantee  existence  of  a  robust  controller  for  the 

A  A  A 

modified  system,  Davison's  results  in  terms  of  A,  B  and  E 
can  be  restated  as  follows: 

Theorem  1 ' 

If  the  disturbance  w  is  measurable  then  there  exists  a 

feedforward  controller  for  the  augmented  system  (29,33)  so 

that  e(t)  =  (y^  n  -  y)  0  as  t -*  for  all  constant 

ol  r 

disturbances  w  e  R  ,  and  constant  reference  inputs  R 

i f  and  only  if 

rank  G'  (0)  =  rank  C(sl-A)**1  b!  =  r 

I  s=o 

Theorem  2' 

There  exists  a  robust  controller  for  the  augmented 
system  (29,33)  such  that  e(t)-»  0  as  t -•  «>  for  all  constant 

j 

disturbances  w  6  R  and  all  constant  reference  inputs 
y e  ^  anc*  t^ie  c^osec*  looP  system  is  asymptotically 
stable  if  and  only  if 

rank  G'  (0)  =  rank  C(sl-A)  '  B  j^_Q  =  r 
Theorems  and  2'  are  elementary  extensions  of 
Davison's  results  [1]. 

Having  ascertained  the  existence  of  a  robust  feedback, 
feedforward  controller,  the  final  controller  is  given  by: 


' 
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u  =  t  T !  ^  V  v-c  f  “  [  T  i'  ]  ;  D  [  w  -  [l\  ]  '  £  J  {y-y^f  )dV  (34) 

where  T  j  =  G'(0)  and  D|  =  C(sl-Af‘  E  a'nd  £>  0. 

It  is  interesting  to  note  that  block  diagram 
manipulation  of  this  scheme  shows  that  the  overall  'robust7 
control  scheme  is  in  fact  a  multivariable  PI  scheme  as  can 
be  seen  from  figure  6.  Notice  also  the  different  amount  of 
proportional  action  taken  on  the  y^ef  and  y  signals.  From  a 
frequency  domain  viewpoint  the  controller  matrix  [T('  ]  1  is 
simply  a  steady-state  decoupling  matrix. 


2.5  Evaluation  of  Robust  Controller  on  a  Simulation  Program 

2.5.1  Introduction 

In  this  chapter  we  will  describe  the  design  and 
implementation  of  a  robust  feedback  feedforward  controller 
on  a  nonlinear  distillation  column  simulation  program. 

First  we  will  have  a  look  at  the  program  itself  and  the 
way  the  slave  loop  control  was  implemented. 

Then  we  will  design  a  robust  control ler  for  the 
simulated  plant,  and  finally  we  will  tune  the  robust 
controller  and  compare  its  performance  to  a  conventional  PID 
control ler . 

2.5.2  The  Distillation  Column  Simulation  Program 

The  simulation  package  that  was  used  in  this  study  is 
an  extension  of  the  original  programs  that  were  developed  by 
Svrcek  [17]  to  model  the  department's  pilot  distillation 
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Figure  6 
slave  loops 


Robust  controller  with  combined  master  and 
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co 1 umn . 

The  model  is  nonlinear  and  based  on  material  and  energy 
balances.  Because  the  model  presents  an  adequate 
representat ion  of  binary  distillation  column  behaviour,  it 
is  the  most  useful  tool  for  evaluating  controller  behaviour 
short  of  actual  plant  implementation.  Due  to  changes  in  the 
equipment  and  operating  conditions  of  the  pilot  scale 
distillation  column  (described  in  part  4  of  this  thesis), 
which  are  not  yet  reflected  by  the  simulation  programs,  we 
do  not  claim  that  the  model  accurately  presents  the  dynamics 
of  the  pilot  plant  distillation  column.  Therefore  the 
simulations  and  experimental  evaluation  should  be  seen  as 
two  separate  entities. 

In  figure  7  we  see  the  proposed  implementation  of  the 
robust  controller.  Note  that  the  proportional  controller  in 
the  slave  loop  is  situated  in  the  feedback  loop  and  not  in 
the  forward  path  as  originally  proposed  in  section  2.4.  The 
reasons  for  this  are: 

1.  On  a  real  plant  one  could  have  a  conventional 
controller,  traditionally  located  in  the  forward  path, 
which  provides  analog  backup  in  case  of  a  computer 
failure.  The  inputs  to  controllers  in  the  forward  paths 
would  be  set  points.  (In  our  simulation  program  we  omit 
analog  controllers  for  backup.) 

2.  The  objective  of  the  simulations  was  not  only  controller 
evaluation,  but  also  familiarization  with  the 
distillation  column  dynamics.  In  the  open- loop  case  the 
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Figure  7.  Robust  controller  for  simulation  program 

control  inputs  for  the  distillation  column  are 
distillate  flow  rate  (top  product)  and  reboiler  steam 
flow  rate.  In  the  proportionally  controlled  closed  loop 
case  with  conventional  controllers  in  the  forward  paths, 
the  inputs  are  set  points  for  top  product  composition 
and  bottom  product  composition.  With  the  controllers  in 
the  feedback  loops,  the  inputs  will  remain  in  the  same 


i  *< 


y 


32 


engineering  units  (distillate  flow  rate  and  steam  flow 
rate,  as  in  the  open  loop  case).  We  felt  that  this 
allows  us  to  better  demonstrate  the  linearizing  effect 
of  proportional  control  and  also  the  responses  of  top 
and  bottom  composition  to  distillate  flow  rate  and  steam 
flow  rate  rather  than  set  points  gave  us  a  better 
insight  in  the  column  dynamics,  especially  the 
interaction  patterns. 

Note  that  in  this  figure  Tj  will  be  different  from  t[  in 
figure  5. 

In  the  simulation  programs,  the  control  loops  are: 
top  loop:  input  =  distillate  flow  in  grams/sec. 

output  =  top  product  composition  in 
weight  fraction  methanol 

bottom  loop:  input  =  reboiler  steam  flow  rate  in 

grams/ sec . 

output  =  bottom  product  composition  in  . 
weight  fraction  methanol 
The  steady  state  operating  conditions  are: 

Feed  :  flow  rate  18.055  g/sec,  50%  methanol 

Distillate:  flow  rate  9.025  g/sec,  95.986  %  methanol 
Bottom  :  flow  rate  9.030  g/sec,  4.280  %  methanol 
:  flow  rate  14.577  g/sec 


Steam 


hr  .  i  *h  I  f  rb<  O  fT*  t  •ip'^ 
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2.5.3  Controller  Design 

As  was  outlined  in  section  2.2,  the  robust  control ler 
was  designed  for  linear  systems.  On  certain  Kinds  of 
nonlinear  systems  we  have  to  use  proportional  feedback 
control  first  to  linearize  the  system  behaviour  (see  section 
2.4).  We  will  try  to  design  a  controller  for  the 
distillation  column  model  and  show  the  differences  between  a 
plant  with  and  without  a  proportional  slave  loop. 

The  experiments  that  are  required  for  the  controller 
design  were  first  performed  on  the  open- loop  plant  model: 

Disturbance  Top  composition  Bottom  composition 


+5% 

feed  = 
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For  positive  disturbances, 


'-9. 80  -9.80 

.  89.7  -55.9 


det  [ T ‘ ]  =  0.00070 
For  negative  disturbances, 

f-  1 0 . 2  -10.2 

[T|]  = 

1.86.4  -16.3 


det  [ T | ]  =  0.00096 

Although  the  determinant  of  the  t| -matrix  for  negative 
changes  is  37%  larger  than  for  positive  changes,  the 
controller  matrices  seem  to  be  not  too  different.  There  is 
however  a  significant  difference  in  element  (2,2)  of  the 
[ T ' ]  1  mat r i ces . 

The  problem  we  indicated  in  section  2.3,  that  of  sign 
changes  in  the  controllers,  will  not  arise  in  this  control 
setup  because  the  [t|]  matrices  are  positive  definite: 


a  matrix 


wi  th 


. alz  positive 


-  an  ~a2i 

will  always  have  a  positive  determinant,  regardless  of  the 
magnitude  of  the  entries.  To  change  sign,  one  of  the  signs 
of  the  entries  would  have  to  change  which  means  that  the 
direction  of  one  of  the  responses  would  have  to  change.  This 
will  not  happen  for  normal  operating  conditions. 

The  only  problem  we  encounter  here  is  a  large  uncertainty  in 
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the  magnitude  of  the  entries,  but  not  in  the  signs.  In 
another  part  of  this  thesis  we  will  see  examples  of  more 
problematic  plant  behaviour. 

To  avoid  the  problems  described  above,  we  apply 
proportional  feedback  control  first  (see  fig.  8) 

The  gains  for  top  and  bottom  loops  (k,  and  kt )  were 
chosen  somewhat  arbitrarily  as  95.0  and  80.0.  Some  runs  for 
feed  flow  disturbances  showed  that  these  gains  gave  fast  and 
smooth  responses  without  oscillations  (see  appendix  B  plots 
3  and  4 ) . 

The  controller  design  experiments  on  the  proportionally 
controlled  plant  yield: 

For  positive  changes: 


-i 


’'-  105  -9.72' 

.87.2  -110. 


0.00366 

0.01100 


For  negative  changes: 


-105 

-9.91' 

0.00388’ 

D!  = 

112 

-104  . 

0.00975, 

As  we  can  see,  the  differences  between  positive  and 
negative  control ler  matrices  are  fairly  small,  especially  in 
the  diagonal  terms,  so  the  objective  of  "plant 
linearization"  has  been  achieved  to  our  satisfaction. 


The  control ler  we  will  use  is: 

A 


u  =  [T|  f  {y 


-  D>  -  £  (y  -  y,tf  )dt'  -  Kpy 


(35) 
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composition 


Figure  8 


where : 


Plant  with  slave  control  loop 


1 05  -9.80' 

0.00377' 

d!  = 

99.5  -107 

J 

0.01040 

the  average  of  positive  and  negative  experimental  results, 

95  O' 

0  80 

m 

This  controller  was  implemented  in  a  discretized  version 
with  sample  time  128  secs. 

2.5.4  Robust  Feedback -Feedforward  Control 

The  next  stage  after  design  and  implementation  of  the 
controller  is  the  tuning,  done  by  manipulating  epsilon  in 
equation  (35)  .  The  tuning  was  done  for  a  -20%  feed  flow 
disturbance  because  the  primary  objective  is  regulatory 
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control  of  the  distillation  column. 

There  are  well  tuned  conventional  multi  loop  PI  and  PID 
controllers  available  that  are  tuned  for  the  same 
di sturbance . 

The  tuning  and  performance  comparison  was  based  on  the 
calculation  of  the  Integral  Absolute  Error  (IAE)  for  192 
minutes  (96  samples). 

For  robust  control  without  feedback  ( D[  =  0)  the  tuning 
results  are  listed  in  table  1. 

The  best  results  were  achieved  for  epsilon  =  0.0008 
If  we  include  feed  forward  action  in  the  robust 
controller,  the  tuning  results  are  as  outlined  in  table  2. 

We  see  that  in  this  case  the  control  is  best  if  the 
robust  controller  has  feed  forward  control  only.  Since  the 
distillation  column  model  is  completely  noise  free  and  the 
parameters  are  time  invariant,  feed  forward  control  is 
ideally  suited  to  handle  load  disturbances.  Obviously  we 
have  to  compromise  to  get  realistic  results.  We  chose 
epsilon  =  0.0003  to  continue  the  simulations  for  the  robust 
feedback  feed  forward  controller. 

2.5.5  Comparison  to  PI  controllers 

To  evaluate  the  performance  of  the  robust  controller  we 
use  well  tuned  conventional  multi  loop  PI  controllers  on  the 
same  plant  model  for  the  same  operating  conditions. 


' 
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Table  1 . 


Tuning  of  robust  feedback  controller 


IAE 


eps i Ion 

top 

bot  tom 

0.0006 

23.3 

62.8 

0.0008 

20.9 

48.7 

0.0010 

23.6 

49.8 

Table  2 . 


Tuning  of  robust  feedback- feedforward  controller 


IAE 


epsi Ion 

top 

bot  tom 

0.0008 

15.  1 

32.6 

0.0004 

10.2 

22.9 

0.0003 

9.7 

22.2 

0.0002 

9.2 

21  .6 

0.0000 

7.4 

17.2 

The  PI  controllers  were  tuned  for  a  -20%  feed  flow 
disturbance,  the  controller  constants  are: 

top  loop:  Kp=  95,  integral  time  Tx  =  1000  secs, 

bottom  loop:  Kp=  80,  integral  time  Tx  =  1450  secs. 

Table  3  gives  a  comparison  between  robust  feedback  and  PI 
control . 

Appendix  B  plots  5  to  10  show  some  responses  for  robust 
feedback  and  PI  control. 


To  evaluate  the  robust  feedforward  feedback  control ler  we 
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Table  3. 

Robust  feedback  control  vs.  PI  control 

robust  PI 


disturbance 

top 

bottom 

top 

bot  tom 

+20%  feed  flow 

17.  1 

51  .4 

19.0 

57.5 

-20%  feed  flow 

20.9 

48.7 

21  .5 

48.9 

+1%  top  composition 

4.4 

6.8 

6.2 

25.5 

-1%  top  composition 

9.3 

12.5 

6 . 6 

13.0 

+1%  bottom  composition 

5.4 

14.3 

3.7 

10.6 

-1%  bottom  composition 

4.0 

10.8 

2.9 

8.6 

61.1 

144.5 

59.9 

164.  1 

205.6  224.0 

used  a  P I -p 1  us- feedforward  controller  for  comparison.  In 
this  controller  we  used  the  same  type  of  feedforward  action 
as  in  the  robust  controller,  i.e.  we  added  to  the  multi  loop 
PI  controllers  the  term  - [ Tj  T’dJw,  where  w  is  the  feed  flow 
rate  and  D  and  [Tj  j”1  are  the  same  as  in  equation  (35)  .  The 
performance  comparison  is  given  in  table  4. 

Appendix  B  plots  11  and  12  show  responses  for  -20%  feed  flow 
di s turbances . 

In  the  robust  feedback  case  we  find  that  the  robust 
controller  is  slightly  better  than  the  PI  controller  except 
for  set  point  changes  in  the  bottom  composition.  The  overall 
control  performance  however  is  slightly  better  for  the 
robust  controller. 
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Table  4. 

Robust  feedback- feedforward  vs.  Pl+feedforward 


robust 

PI 

disturbance 

top 

bot  tom 

top 

bottom 

+20%  feed  flow 

6.5 

21  .7 

7.2 

22.5 

-20%  feed  flow 

3.8 

11.7 

13.1 

29.7 

+1%  top  composition 

3.8 

11.7 

6.2 

25.5 

-1%  top  composition 

6.3 

11.2 

6.6 

13.0 

+1%  bottom  composition 

3.8 

11.4 

3.7 

10.6 

-1%  bottom  composition 

3.6 

9.8 

2.9 

8.6 

33.7  88.0 

39.7 

109.9 

121  .7 

149 

.6 

If  we  add  feedforward  action  to  the  controllers,  as 
shown  in  table  4,  the  gap  widens.  The  robust 
feedback- feedforward  controller  performs  better  than  the 
P I  -  feedforward  controller  for  all  load  and  set  point  changes 
except  set  point  changes  for  the  bottom  composition,  where 
the  difference  is  minimal. 

If  we  retune  the  PI  controller  after  adding  feedforward 
control,  we  find  that  the  controller  performance  for  a  load 
change  improves  with  increasing  integral  times  (i.e. 
decreasing  integral  action).  In  fact  the  absence  of  noise 
and  parameter  drift  in  the  plant  model  allows  us  to  drop  the 
integral  action  completely  if  we  have  a  good  feedforward 
controller.  But  then  the  PI  controller  performance  for  set 
point  changes  deteriorates  sharply.  We  have  not  found  any 
one  PI  controller  with  an  overall  performance  surpassing  the 
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robust  controller  performance. 

2.5.6  Addition  of  Derivative  Action 

It  might  be  possible  to  improve  the  controller 
performance  by  adding  derivative  action  to  tthe  controller. 
Now  the  control  algorithm  (35)  will  be  modified  to: 


D',w  -  £  (y  -  y^ej>  )dt #  } 


O 


(36) 


and  =  (k,  Tjt|  )  / 1  SEC ,  Tj  is  the  derivative  time  for  the 
top  loop 

We  basically  used  the  same  robust  controller  as  in  section 
2.5.5,  but  we  added  the  derivative  action  and  tuned  the 
derivative  times  for  a  -20%  feed  flow  disturbance.  For  this 
case  the  optimum  derivative  times  were  110  secs,  (top  loop) 
and  180  secs,  (bottom  loop).  For  comparison  we  used 
multi  loop  PID-plus-feedforward  controllers,  tuned  for  a  -20% 
feed  flow  disturbance: 

top  loop:  Kp=  95  Tx  =  2500  secs  =  80  secs, 

bottom  loop:  Kp=  80  TL  =  3000  secs  Tj  =120  secs. 

See  table  5  for  the  performance  comparison. 

Appendix  B  plots  13  to  18  show  load  disturbances  and  set 
point  changes  for  PD-robust  feedback  feedforward  control  and 
for  PID-plus-feedforward  control. 

Again  we  see  that,  as  before,  the  overall  performance  of  the 
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Table  5. 

PD-robust  +  FF  vs.  PID+FF 

robust  PID 


disturbance 

top 

bottom 

top 

bot  tom 

+20%  feed  flow 

5.9 

19.0 

6.  1 

18.9 

-20%  feed  flow 

7.  1 

17.4 

7.9 

18.1 

+1%  top  composition 

4.4 

11.8 

5 . 0 

45.2 

-1%  top  composition 

5.5 

10.2 

3.6 

22.5 

+1%  bottom  composition 

2.7 

10.8 

2.9 

9.  1 

-1%  bottom  composition 

2.7 

9.7 

2.5 

10.2 

28.3 

78.9 

28.0 

124.0 

107 

.2 

152 

.0 

robust  control ler  is  better  than  the  conventional  PID 
control ler . 

2.5.7  Conclusions 

The  previous  sections  proved  that  the  thesis  "Robust 
control  is  at  least  as  good  as  or  better  than  conventional 
P I ( D )  control"  is  justified,  even  for  a  nonlinear  plant. 

For  the  more  important  load  changes,  which  would  occur 
in  regulatory  controlled  situations,  the  robust  controller 
surpassed  the  conventional  controller  in  all  simulations. 

For  set  point  changes  we  saw  that  for  top  composition 
changes  generally  the  robust  controller  was  better  whereas 
for  bottom  composition  changes  the  P I ( D )  controller  was 
slightly  better. 


- 
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The  fact  that  the  simulation  programs  are  completely 
free  of  noise  and  parameter  drift  posed  a  problem  in  the 
tuning.  A  well  designed  feedforward  controller  will  bring 
the  output  variables  so  close  to  their  set  points  that 
integral  action  is  hardly  necessary  to  reduce  the  offset, 
or,  the  best  integral  action  is  no  integral  action. 

Rather  than  spending  much  effort  and  time  on  building 
noise  and  parameter  drift  into  the  simulation  programs  we 
decided  to  use  the  results  as  they  are  because: 

1.  we  rather  spend  some  extra  time  on  the  more  realistic 
pilot  plant  evaluation  as  outlined  in  chapter  4. 

2.  Both  robust  and  P I ( D )  controllers  suffered  from  the  same 
problem  so  the  comparison  is  still  valid. 

Not  reflected  in  the  simulation  results  is  the  tuning  effort 
which  was  considerably  less  for  the  robust  controller  with 
only  one  tuning  constant  (epsilon)  vs.  the  PI  controller 
with  4  tuning  constants. 

The  results  as  presented  in  this  chapter  justify  and 
strongly  encourage  a  more  realistic  pilot  plant  evaluation 
of  the  robust  controller. 
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3.  Development  of  a  Control  System 


3.1  Introduction 

In  order  to  apply  the  robust  controller  as  described  in 
the  previous  chapter  to  a  pilot  plant  distillation  column  we 
need  to  develop  a  control  system,  i.e.  putting  the  hardware 
together  and  developing  the  software.  Since  small  and 
powerful  1  computers  have  become  available  there  has  been  a 
tendency  to  use  these  mini  or  micro  computers  in  a  network 
to  increase  flexibility  and  reliability  [19,20,22,23]. 

This  chapter  describes  the  requirements  for  a  process 
control  system  of  the  type  involved  in  this  study.  Both 
hardware  and  software  requirements  of  the  computer  control 
system  are  considered. 

In  section  3.2  the  hardware  configuration  of  the 
control  system  is  described,  together  with  how  it  ties  in 
with  the  existing  instrumentation  and  the  distillation 
column.  This  is  followed  by  a  description  of  the  logical 
setup  of  the  complete  controlled  system,  i.e.  which  tasks 
are  situated  in  which  computer  and  why.  The  development  of 
the  mi crocomputer  software  is  described  next. 

The  software  control  system  has  been  designed  to  be 
flexible  and  modular.  This  aspect  of  the  system  design  is 
also  described  in  this  chapter,  i.e.  the  amount  of  effort 
involved  in  making  changes  in  the  programs.  The  chapter  is 
concluded  with  an  evaluation  of  the  newly  implemented  total 

control  setup. 
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3.2  Requirements  for  the  Control  System 

It  is  important  to  consider  the  main  requirements  for  a 

control  system  software  package  for  a  University  research 

environment  as  opposed  to  an  industrial  system  [18-20]. 

The  main  requirements  are: 

perform  simple  conventional  control  (PID); 
provide  sufficient  data  acquisition  capabilities: 
provide  an  operator -process  interface  via  a  CRT  console 
easily  changeable  control  subroutines; 
no  limitations  (except  program  size)  on  the  complexity 
of  control  algorithms  that  can  be  implemented,  i.e.  any 
algorithm  that  is  Fortran  programmable  and  does  not 
require  excessive  computer  memory  can  be  used; 
easy  on-line  changeability  of  I /O-parameters  in  case  of 
a  system  reconf igurat ion . 

The  differences  with  an  industrial  system  are  mainly: 

1.  An  industrial  system  is  generally  limited  to  a  number  of 
fixed  control  algorithms. 

2.  Industrial  systems  place  heavy  emphasis  on  alarm  display 
and  acknow 1 edgemen  t . 

3.  In  an  industrial  system  the  operator  is  not  required  to 
make  on-line  changes  in  I /0-conf i gurat i on . 

4.  Operator  interfacing  is  done  through  functional  push 
buttons  on  a  specially  designed  console  rather  than  a 
regular  teletype  keyboard  and  a  terminal. 

The  requirement  of  on-line  changeability  of  I /O-parameter s 

also  assumes  that  all  instrumented  measurement  signals  are 
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physically  connected  to  the  computer's  I/O  facilities  so 
that  reading  those  signals  requires  only  software  changes. 

3.3  Hardware  Configuration  of  Control  System 

To  get  an  understanding  of  the  system  configuration  we 
need  to  consider  the  instrumentation  of  the  distillation 
co 1 umn . 

All  flows  are  determined  by  measuring  differential 
pressures  across  orifice  plates.  The  differential  pressure 
signals  are  converted  to  3-15  psi  air  signals  by  d/p-cells 
(differential  pressure  cells).  The  3-15  psi  air  signals  are 
normally  routed  to  both  pneumatic  analog  controllers  and  P/I 
(air  pressure  to  electric  current)  converters.  The  currents 
from  the  P/I  converters  are  converted  to  voltages  which  are 
hooked  up  to  the  analog  inputs  of  the  LSI-11. 

All  temperatures  are  measured  by  (type  J) 
thermocouples.  The  electrical  signals  from  the  thermocouples 
are  connected  to  both  a  Honeywell  temperature  recorder  and  a 
low  level  analog  input  module  in  the  LSI-11. 

Pressures  and  liquid  levels  are  measured  with 
d/p-cells.  The  3-15  psi  air  signals  from  the  d/p-cells  are 
connected  to  P/I  converters.  The  currents  from  the  P/I 
converters  are,  after  conversion  to  voltages,  connected  to 
the  LSI-11  analog  inputs. 

The  top  product  composition  is  measured  with  a  gas 
chromatograph  (GC)  which  automatically  samples  once  every 
three  minutes.  The  GC-analysis  is  being  done  by  an  HP  1000 
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computer.  This  computer  produces  a  GC-report  which  is  sent 
to  the  LSI-11  (through  an  RS-232  line). 

All  controlled  variables  on  the  distillation  column  are 
flow  rates.  These  flow  rates  are  controlled  by  conventional 
pneumatic  control  valves  which  get  their  input  signals  from 
analog  controllers. 

The  set  points  for  the  analog  controllers  are  pneumatic 
signals  from  either  manual  valves  in  case  of  a  local  flow  or 
level  control  loop,  or  from  I/P  (electric  current  to  air 
pressure)  converters  in  case  of  computer  controlled  loops. 
The  I/P  converters  get  their  inputs  (electrical  currents 
10-50  mA )  from  Current  Output  Stations  (COS)  which  in  turn 
get  their  inputs  (0-10  Volts)  from  the  analog  outputs  of  the 
LSI-11.  See  figure  9  for  a  schematic  diagram  with  one 
control  loop. 

The  LSI- 11  is  a  16  bit  word  length  microcomputer  with 
28k  words  of  read/write  memory  [21].  This  mi crocomputer 
forms  the  center  of  a  fairly  powerful  system  which  also 
contains  10-devices  (interface  modules  for  high  level  analog 
in,  low  level  analog  in,  digital  in,  digital  out,  analog  out 
and  serial  line  interface  modules)  and  mass  storage  devices 
(floppy  disk  system). 

Compared  to  a  conventionally  controlled  plant,  the 
LSI- 11  microcomputer  replaces  analog  controllers  and  chart 
recorders.  It  also  provides  new  possibilities  that  do  not 
exist  in  a  conventionally  instrumented  plant  or  a  loop,  such 
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Figure  9.  Distillation  column  with  instrumentation  for 
top  loop 

as  composition  control  on  the  bottom  loop  since  the 
mi crocomputer  can  handle  an  ASCII  GC-report  (see  next 
section),  that  a  pneumatic  controller  cannot.  Another  new 
possibility  is  the  use  of  complex  control  algorithms  which 
are  very  difficult  to  realize  or  implement  using  analog 
control lers . 

The  total  cost  of  the  computer  system  is  approximately 
C$  12,000  (1979).  This  number  includes  CPU,  memory,  a  dual 
floppy  disk  drive  and  all  process  and  terminal  interfaces. 


; 
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3.4  Overall  Data  Acquisition  and  Control  System 

The  overall  configuration  of  the  microcomputer  control 
system  is  shown  in  figure  10.  This  figure  is  also  useful  in 
describing  the  various  measurement  systems  and  links  between 
the  mi crocomputer  and  the  HP- 1000. 

The  distillation  column  with  the  LSI-11  system  can  be 
considered  as  a  totally  independent  plant.  The  only 
"external"  device  needed  is  the  GC  computer  for  the  bottom 
product  composition  measurement.  Any  failure  in  any  of  the 
other  multi-user  computer  systems  will  not  affect  the 
distillation  column. 

The  LSI- 11  takes  care  of  all  process  input,  control  and 
process  output,  operator  interfacing  and  historical  data 
storage.  Except  for  the  bottom  product  composition  all 
process  inputs  are  read  directly  from  the  distillation 
column  instrumentation.  All  output  signals  are  sent  directly 
to  the  distillation  column  instruments. 

The  bottom  product  composition  is  determined  from  a  gas 
chromatogram  by  the  GC-computer.  This  computer  produces  a 
GC-report  which  is  intended  to  be  printed  on  a  normal 
terminal.  For  our  control  applications  the  GC-report  is 
being  sent  to  a  serial  line  interface  module  in  the  LSI-11. 
Then  the  report  will  be  analyzed  to  find  the  bottom  product 
compos i t i on . 

The  data  acquisition  is  set  up  to  send  one  line  of  data 
per  sample  interval  to  a  floppy  disk  file.  The  line  of  data 
consists  of  the  sample  number,  time,  relevant  process 
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Figure  10.  Schematic  diagram  of  the  distributed  network 
of  mini  and  micro  computers  and  the 
interface  with  the  computer  controlled 
distillation  col umn 
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readings  and  controller  outputs.  It  is  possible  to  store  all 
data  for  two  weeks  of  uninterrupted  column  operation  on  one 
floppy  disk  (approximately  900  blocks  of  128  words  each). 

There  exist  data  transmission  facilities  between  the 
LSI-11  and  the  higher  level  HP  1000  minicomputers.  This  line 
can  be  used  to  transmit  data  to  the  HP  1000  for  plotting  and 
report  generating.  Future  applications  of  supervisory  and 
optimizing  control  at  a  higher  level  are  also  possible. 

There  is  also  a  data  transmission  link  between  the 
department's  HP  system  and  the  University's  Amdahl 
maxicomputer  which  allows  even  more  sophisticated  plotting 
and  report  generation. 

In  previous  studies  IBM  1800  and  HP  1000  computers  were 
used  for  control  ,  in  the  latter  case  the  process  10  was 
handled  by  an  HP-Macs  system.  Because  of  the  fact  that  the 
LSI- 11  is  a  dedicated  system  whereas  the  previous  ones  were 
multi-user  systems,  the  reliability  has  significantly 
increased  since  there  is  no  danger  anymore  of  other  users 
bringing  the  computer  down. 

For  safety  reasons  there  is  still  a  "crash  program" 
running  in  the  HP  1000  system  which  takes  readings  of  all 
relevant  process  variables,  parallel  to  the  LSI-11  (not 
shown  in  figure  10).  If  variables  exceed  safe  operating 
limits,  the  crash  program  automatically  shuts  down  the 
distillation  process  by  tripping  the  main  power  and  air 
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In  the  previous  paragraphs  we  saw  that  the  "lower  level 
tasks"  (control  and  data  acquisition),  which  require  a  high 
degree  of  reliability,  are  located  in  a  totally  dedicated 
microcomputer ,  which  also  takes  care  of  presentation  of  data 
to  the  operator,  "Higher  level"  off-line  tasks  such  as 
appendix  B  plotting  and  report  generation  that  do  not 
require  a  high  integrity  system  are  located  in  a  more 
sophisticated  multi-user  system  that  offers  a  wider  array  of 
user  services.  In  this  set  up  we  tried  to  optimize  the  use 
of  computer  facilies  to  obtain  integrity,  flexibility  and 
sophisticated  user  facilities,  the  main  reasons  for  the 
current  trend  towards  distributed  computer  systems  [22,23], 


3.5  Development  of  the  Control  Software 
3.5.1  Introduction 

To  meet  all  the  requirements  for  the  control  system, 
the  microcomputer  has  to  perform  two  different  types  of 
tasks : 

a.  tasks  directly  related  to  process  control  and  data 
acquisition  such  as  input,  control  calculation  and 
send  the  valve  signals; 

b.  tasks  related  to  operator  interfacing,  i.e.  the 
presentation  of  data  to  the  operator  and  the 
changing  of  parameters  by  the  operator . 

The  major  differences  from  a  software  point  of  view  are  that 
the  process  control  tasks  are  time  critical,  they  should  be 
pgp formed  every  sample  time  at  the  right  moment  and  in  the 
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proper  sequence. 

On  the  other  hand,  the  operator  should  be  able  to  enter 
changes  in  the  parameters  whenever  he  likes  and  also  list 
data  whenever  he  likes. 

Also,  the  process  control  and  data  acquisition  cycle  should 
have  a  high  degree  of  integrity,  maintaining  a  controlled 
situation  in  a  plant  is  more  important  than  data 
presentation  to  the  operator.  Operator  mistakes  should 
affect  the  control  performance  as  little  as  possible. 

To  achieve  the  desired  control  integrity  and  operator 
interface  independence,  several  of  the  advanced  features  the 
LSI-11  and  its  operating  system  (RT-11)  are  described 
[24-26] . 

The  software  was  written  in  two  parts,  called  a 
foreground  job  and  a  background  job.  (RT-11  feature).  See 
also  figure  1 1 . 

The  high  priority  foreground  job  consists  of  a  mainline 
foreground  program  and  several  subroutine  libraries.  Its 
tasks  are  to: 

perform  the  actual  control,  i.e.  read  inputs,  calculate 
control  and  send  outputs; 

take  care  of  data  acquisition  on  floppy  disks; 
perform  the  program  timing,  i.e.  keep  track  of  the  time 
and  schedule  the  control  cycle  every  sample  time. 

The  lower  priority  background  job  consists  of  a 
mainline  background  program  and  several  subroutine 
libraries.  Its  tasks  are  to  take  care  of  the  operator  input 
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F igure  1 1 . 


Overall  flow  diagram  for  foreground  and 
background  programs 
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and  output,  that  is  to: 

print  relevant  process  data  every  sample  time  and/or  on 
demand ; 

interpret  and  implement  operator  commands. 

Although  there  is  communication  between  the  background  and 
foreground  programs,  it  is  not  necessary  for  the  programs  to 
run  simultaneously.  The  setup  was  chosen  this  way  to  enable 
the  operator  to  issue  commands  and  interpret  IO-without 
interfering  with  the  process  control.  An  operator  I/O 
mistake  that  would  cause  the  program  to  abort,  will  in  this 
case  only  abort  the  background  program  without  affecting  the 
control  performance. 

We  will  first  look  at  the  way  data  is  handled  in  the 
programs  and  the  use  of  common  areas.  This  will  be  followed 
by  a  detailed  description  of  the  foreground  and  background 
programs.  Data  acquisition  on  floppy  disks  will  be  discussed 
in  the  next  section  and  we  will  give  an  overview  of  the 
operator  commands.  The  last  section  will  be  a  guideline  for 
future  users  who  want  to  tailor  the  system  to  their  own 
needs . 

3.5.2  Buffer  Usage  in  Foreground  and  Background  Programs 

Most  of  the  data  transmission  between  the  programs  and 
subroutines  is  being  handled  by  labelled  COMMON  areas.  The 
reasons  for  using  COMMON  instead  of  subroutine  arguments 

are : 

1.  saving  of  memory  space; 
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2.  all  parameters  and  data  for  transmission  have  to  be 

assembled  in  an  integer  buffer.  To  put  real  values  in  an 
integer  buffer  we  use  EQUIVALENCE  statements.  Since 
EQUIVALENCEd  variables  may  not  appear  as  subroutine 
arguments,  one  has  to  use  COMMON. 

There  are  four  COMMON  areas: 

BLK1 :  The  I/O-table..  This  is  the  table  containing  all 
Input-Output  vectors. 

Dimension  I0TAB(8,12) 

I /O-vectors:  Each  of  the  12  I/O-vectors  in  the  I/O-table 
is  8  words  long.  The  layout  of  a  vector  is: 


Mn£M  OMtc 

VALUE 

INFO 

S  P  A  A/ 

2 

MNEMONIC:  (words  1  and  2)  4  characters,  e.g.  TOPC  for 
TOP  product  Composition.  An  unused  vector  starts  with 
the  characters  UN.  Any  vector  starting  with  UN  will  be 
ignored  by  input  and  output  routines. 

VALUE:  (word  3)  an  integer  number  representing 
millivolts.  For  an  output  vector  the  value  would  be 
between  0-10000  (mVolt).  For  an  input  the  value  would  be 
between  1000-5000  (mVolt).  So  the  resolution  for  an 
input  reading  is  1/4000. 

INFO:  (word  4)  Two  bytes  containing  information  about 
the  type  of  I/O-vector.  The  layout  is  as  follows: 


CWA/V/V6X  NUM9b* 


AS  SHo^j/u  Sftow  <Jr 


lb  IS  H/  fi  ll  II  lo  9  S  1  t  5  L,  3  2  t 
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Low  byte:  bit  1=1  for  input,  0  for  output 

2  =  1  for  analog,  0  for  digital 
3<54  =  0  for  linear  input 

1  for  square  input 

2  for  square  root  input 

3  for  temperature  input 
5-8  =  not  used 

High  byte:  bits  9-16  =  channel  number  0-256 
SPAN:  (words  5  and  6)  Real  value  giving  the  span  of  the 
signal.  This  is  necessary  for  conversion  to  engineering 
units. 

ZERO:  (words  7  and  8)  Real  value  giving  the  zero  of  the 
signal,  necessary  for  conversion  to  engineering  units. 
example :  the  vector  for  the  steam  flowrate  input: 
MNEMONIC:  STFL 
VALUE  :  1000-5000 
INFO  :  channel  #0,  SQRT  Analog  In 
bit  pattern  0000000000001011 
SPAN  :  22.5002 
ZERO  :  0.0 

See  the  documentation  of  subroutine  ENG  for  an 
explanation  of  SQRT,  SPAN  and  ZERO. 

BLK2 :  The  controller  parameters  table. 

Dimension  CPTAB(40) 

This  is  a  forty  word  array  containing  all  control 
parameters.  These  are  all  the  parameters  that  are 
required  by  the  control  algorithm. 


,  * 
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The  CP- tab le  is  EQUIVALENCEd  to  the  following 
parameters : 


word  i 

oarameter 

size 

mean i ng 

CPTAB (  1  ) 

I  SEC 

1 

sample  interval  in  seconds 

CPTAB ( 2 ) 

ITYPE 

1 

type  of  control 

CPTAB ( 3 ) 

YREF (2) 

4 

set  points  for  top  and  bottom 

compos i t i on 

CPTAB ( 7 ) 

KP  ( 2 ) 

4 

P-  gains  for  top  and  bottom 

loops 

CPTAB ( 1 1 ) 

T  X  ( 2 ) 

4 

integral  times  for  top  and 

bottom  loops 

CPTAB ( 15) 

TD  ( 2 ) 

4 

deriv.  times  for  top  and 

bottom  loops 

CPTAB ( 19) 

Til  (2,2) 

8 

matrix  for  robust  control 

CPTAB ( 27 ) 

D1  (2) 

4 

feed  forward  control  matrix 

CPTAB ( 3 1 ) 

EPS ( 2 ) 

4 

tuning  constants  for  robust 

control 

CPTAB ( 35 ) 

SI AE ( 2  ) 

4 

integral  absolute  errors  for 

both  loops 

CPT  AB ( 39  ) 

not  used 

CPTAB ( 40 ) 

not  used 

BLK3 :  The  program  parameters 

table 

Dimension 

PPTAB( 10) 

This  is  a 

ten  word 

ar  ray 

containing  program  parameters 

for  the  foreground  progr 

am: 

word 

parameter 

size  meaning 

PPTABI 1 ) 

IT  IME ( 2 ) 

2 

clock  time  in  internal  clock 
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format 

PPTAB ( 3 )  NSAMP  1  sample  number 

PPTAB(4)  I SNOR  E ( 2  )  2  waiting  time  for  GC  in  minutes 

and  seconds 

P P T A B ( 6 )  not  used 


PPTAB ( 10)  not  used 
BLK4:  Data  transmission  flags 

FLAGD  =  1  byte  flag  for  data  transmission 
FLAGP  =  1  byte  flag  for  parameter  transmission 
3.5.3  The  Foreground  Program 

After  program  startup,  the  program  first  enters  the 
initialization  phase  where: 

initial  values  for  all  parameters  will  be  read  from  the 
initialization  file  (INIT.DAT); 

the  operator  enters  the  name  of  the  file  to  be  used  for 
data  acquisition. 

Then  the  program  jumps  into  a  continuous  loop: 

1.  If  any  new  parameter  values  are  sent  by  the  background 
program,  those  new  values  are  inserted  in  place  of  the 
old  values,  e.g.  change  from  P  to  PI  control. 

2.  The  GC  report  which  is  being  sent  by  an  HP  1000  computer 
is  read  and  analyzed.  The  result  of  the  analysis  is  a 
number  representing  the  bottom  composition. 

3.  A  timing  routine  waits  for  the  next  sample  time  before 
continuation,  then, 

4.  the  GC  is  triggered  to  get  a  new  sample,  and 


ri  j 

- 


60 


5.  all  specified  analog  input  points  are  read.  This  reading 
gives  us  numbers  for  all  relevant  process  variables. 

6.  The  control  algorithm  calculates  new  settings  for  the 
control  variables,  typically  reflux  flowrate  and  steam 
f lowrate . 

7.  The  specified  Current  Output  Stations  are  addressed  and 
switched  to  the  proper  settings  to  set  the  control 

va 1 ves . 

8.  The  relevant  data  is  written  to  floppy  disk  file 

9.  All  data  and  parameters  are  transmitted  to  the 
background  program 

10.  While  the  GC  is  analyzing  the  next  sample,  the  program 
will  wait  (approx.  2.5  minutes). 

1 1 .  Return  to  1 . 

Figure  12  illustrates  the  flow  diagram  of  this  foreground 
program. 

3.5.4  The  Background  Program 

After  program  startup,  we  first  enter  the 
initialization  phase: 

initial  values  for  all  parameters  will  be  read  from  the 
initialization  file  (INIT.DAT); 

the  "event  file"  (EVENT.DAT)  will  be  opened  and  the  time 
for  the  first  scheduled  event  will  be  read  (see  section 
3. 5. 6. 6  on  event  scheduling. 

Then  the  program  branches  into  a  continuous  loop  and  does 
the  following  tasks: 
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Figure  12.  Flow  diagram  for  foreground  program 
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1.  Check  for  commands  and  scheduled  events  and  services 
them. 

2.  If  there  is  no  new  data  from  the  foreground  program, 
goes  back  to  1 . 

3.  If  there  is  new  data  from  the  foreground  program,  it 
puts  the  data  in  the  tables. 

4.  Prints  relevant  data  on  the  operators  console. 

5 .  Returns  to  1 . 

The  form  of  the  algorithm  is  illustrated  in  figure  13. 

3.5.5  Data  Acquisition  on  Floppy  Disks 

Every  sample  time,  after  the  input-control -output 
algorithms,  the  foreground  program  writes  relevant  data  to 
floppy  disk.  First,  however,  a  disk  file  must  be  defined, 
and  this  can  be  done  in  two  ways: 

1.  Immediately  after  starting  the  foreground  program,  the 
program  will  print  an  asterisk  (*)  in  the  left  margin  on 
the  console.  The  operator  is  expected  to  enter  a 

f i lename : 

First  enter  a  CTRL-F  to  direct  data  to  the  foreground, 
then  enter  the  filename,  e.g.  DY1iDA3132.DAT  (DAta  file, 
opened  on  3  -  13,  2nd  file  of  the  day)  The  operator  can 
of  course  enter  any  name  he  likes. 

Entering  TT:  in  place  of  the  filename  will  direct  the 
data  to  the  console. 

2.  While  the  program  is  running.  By  adding  20  to  the  value 
of  I  TYPE  (see  section  3. 5. 6. 4),  the  operator  can  tell 


* 


' 

m 


63 


Figure  13.  Flow  diagram  for  background  program 
the  foreground  program  to  close  the  old  data  file  and 

open  a  new  one.  ITYPE  will  be  automatically  reset  to  its 

old  value. 

Example:  say  I T Y P E  =  3  and  the  operator  wants  to  close 
the  old  data  file  and  open  a  new  one.  By  doing  so  he 
saves  the  data  in  the  old  file.  The  operator  sets 
ITYPE  to  3+20  =  23.  Just  before  next  sample  moment, 
the  following  will  appear  on  the  console: 

FILE  CLOSED,  ENTER  NEW  FILENAME 
* 

The  operator  replies  CTRL-F  "filename" 

Since  the  foreground  program  hangs  up  while  waiting 
for  the  filename,  it  is  advisable  to  enter  the  new 
filename  immediately  to  prevent  losing  a  sample. 
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Every  sample  time  one  record  of  data  is  being  written 
to  the  disk  file,  this  record  presently  consists  of  11 
integer  numbers: 


word 
1  .  = 
2.  = 

3.  = 

4.  = 

5.  = 

6.  = 

7.  = 

8.  = 

9.  = 

10.  = 
11.= 


sample  number 

time,  minutes  past  midnight 

time,  seconds  past  minutes  past  midnight 

top  product  composition,  a  number  1000  -  5000  (mVolts) 


distillate  flowrate .  1000  -  5000 

reflux  flowrate, .  1000  -  5000 

feed  flowrate .  1000  -  5000 

steam  flowrate .  1000  -  5000 

bottom  product  composition, .  1000  -  5000 

reflux  flowrate  setting, .  0  -  10000 

steam  flowrate  setting, .  0  -  10000 


To  save  different  data  on  file  one  would  have  to  change  the 
'write'  command  in  the  foreground  program. 

3.5.6  Operator  Commands 

Assuming  the  system  is  set  up  properly  and  the  proper 
floppy  disks  are  inserted  in  the  disk  drives,  then  one  can 
proceed  and  run  the  programs  as  follows: 


' 
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3.5.6. 1  Run  the  foreground  program: 

COMMAND:  FRUN  DY 1 : F0R2/N : 2000 
This  command  assumes  that  the  foreground  program  file 
F0R2.REL  is  situated  on  the  floppy  disk  in  drive  1.  The 
additional  space  that  has  to  be  allocated  (in  our  case  1500 
words)  can  be  calculated  from  the  formula  on  page  4-6  of  the 
Advanced  Programmers  Manual  [24]. 

Immediately  after  this  command,  the  program  responds  with  an 
asterisk  (*)  in  the  left  margin.  The  operator  is  now 
expected  to  enter  a  filename  for  data  acquisition.  First 
enter  CTRL-F  (while  pressing  the  control  key,  press  F)  to 
direct  the  filename  to  the  foreground  program. 

Example  (operator  entries  bold) 

* 


CTRL-F . go  to  foreground  program 

DY 1 : DATA . DAT  CR  (Carriage  Return)  filename 
CTRL-B . return  to  background 


Now  the  foreground  program  should  be  running.  The  operator, 
however,  will  not  see  any  listing  on  the  console  until  he 
runs  the  background  program. 

3. 5. 6. 2  Run  the  background  program 
COMMAND:  RUN  DY 1 : BACK3 

This  command  assumes  that  the  background  program  file 
BACK3.SAV  is  situated  on  the  floppy  in  drive  1. 

Note:  One  of  the  things  the  background  program  does  is  to 
take  care  of  the  event  scheduling.  The  program  will  look  for 
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a  file  EVENT.DAT  in  drive  0  and  read  the  first  line  in  this 
file,  i.e.  the  time  for  the  first  scheduled  event.  If  the 
operator  does  not  want  event  scheduling,  the  first  line  in 
the  event  file  has  to  be  00:00,  but  the  event  file  must 
exist.  Look  in  the  next  section  for  details  on  event 
schedu ling. 

Now,  every  sample  time  a  concise  listing  will  appear  on  the 
console.  This  consists  of  two  lines  of  numbers  indicating 
sample  number,  time,  process  variables  and  set  points. 
Before  every  fifth  listing  a  line  of  labels  will  be  printed 
on  the  console: 


TIME 

10:50:31 

TOPC 

95.376 

DIST 

10.318 

REFX 

7.021 

7.203 

FFLO 

18.250 

STFL 

14. 125 
14.100 

BOTC 

5.008 

10:53:31 

95.258 

10. 125 

7.210 

18.105 

14.086 

4.996 

7.312  13.996 

3. 5. 6. 3  The  LIST  command 

There  are  five  printing  options  in  the  background 
program: 

1 .  Concise 

This  consists  of  two  lines  of  numbers  giving  an  overview 
of  the  most  important  process  variables  and  set  points. 

2.  I/O-table 

This  op tion  lists  the  contents  of  the  I/O  table,  i.e.  12 
I/O-vectors.  One  vector  would  look  like: 

TOPC  3000  LIN  A I  1  1  12.8572  89.1673 


■* 

- 

See  the  chapter  on  buffer  usage  (3.5.2)  for  a  detailed 
description  of  I/O-vectors. 
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3.  CP-table 

This  option  lists  the  contents  of  the  CP-table. 

4.  PP-table 

This  option  lists  the  contents  of  the  PP-table 

5 .  All  tab  1 es 

This  option  lists  the  I/O-table,  CP-table  and  PP-table. 
The  operator  can  get  a  listing  on  the  console  by 
entering 

L  return 

and  the  program  will  reply  with  CONCISE(I),  10  TABLE ( 2 )  ,  CP 

TABLE ( 3 ) ,  PP  TABLE ( 4 )  OR  ALL ( 5 ) 

and  the  operator  can  enter  his  choice,  e.g. 

3  return 

This  will  produce  a  listing  of  the  CP-table. 

Instead  of  L  return 

the  operator  could  have  entered  L3  return 

which  would  have  produced  an  immediate  listing  of  the 

CP- table. 

3. 5. 6. 4  The  CHANGE  command 

Instead  of  just  listing  the  table,  the  operator  can 
also  enter  changes  in  the  tables.  To  do  this  the  operator 
has  to  enter  the  numbers  of  the  parameters  he  wants  to 
change  and  the  program  will  then  ask  for  the  new  value  of 
the  indicated  parameter. 
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For  example  the  operator  wants  to  change  the  set  point 
for  the  top  product  composition  from  95%  to  96%.  The  dialog 
is  then: 

C  return 

the  program  replies  with 

CONCISE ( 1 ) f  10-TABLE ( 2 ) ,  CP-TABLE(3),  PP-TABLE(4)  OR  ALL ( 5 ) 
The  set  points  are  controller  parameters,  and  are  to  be 
found  in  the  CP-table.  The  operator  enters 

3  return 

then  the  program  will  list  the  CP-table  and  ask 
ENTER  PARAMETER  NUMBER ( S )  TO  CHANGE: 

In  the  CP-table  listing  the  numbers  are  listed  beside  the 
parameters.  The  set  points  are  parameter  3,  so  the  operator 
enters 

3  return 

and  the  program  will  reply 
NEW  SETPOINTS  YREF : 

Now  the  operator  has  to  enter  both  set  points  (top-bottom): 

96.0,5.0  return 

The  CP-table  will  be  listed  with  the  new  set  points  YREF  and 
the  question 
SATISFIED?  (Y/N) 

Entering  Y  return  will  cause  the  new  parameters  to  be  sent 
to  the  foreground  program  for  implementation. 

Entering  N  return  will  cause  a  reply  by  the  program 
CONCISE(I),  I0-TABLE(2),  CP-TABLE(3),  PP-TABLE(4)  OR  ALL(5) 
then  the  operator  can  enter  the  table  number  and  make 
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additional  changes  on  undo  the  first  change. 

If  the  operator  wants  to  enter  changes  in  the 
I/O-table,  he  will  have  to  enter  two  parameter  numbers, 
indicating  the  position  of  the  parameters  to  be  changed. 
The  listing  of  the  I/O-table  looks  like: 


MNEM. 

VALUE 

INFO 

CH 

SPAN 

ZERO 

1 

TOPC 

3000 

LIN 

A I 

1 1 

12.857 

89. 167 

2 

DIST 

3000 

SQRT 

AI 

8 

15.205 

0.0 

3 

REFX 

1461 

SQRT 

AI 

1 

22.975 

0.0 

4 

FFLO 

3000 

SQRT 

AI 

2 

25.777 

0.0 

5 

STFL 

2400 

SQRT 

AI 

0 

22.500 

0.0 

6 

BOTC 

1625 

LIN 

DI 

0 

32.000 

0.0 

7 

U1SP 

1152 

SQRT 

AO 

2 

22.975 

0.0 

8 

U2SP 

3494 

SQRT 

A0 

3 

22.500 

0.0 

9 

FEED 

5013 

SQRT 

AO 

1 

25.777 

0.0 

10 

UNU2 

0 

LIN 

DO 

0 

100.0 

0.0 

1  1 

UNU3 

0 

LIN 

DO 

0 

100.0 

0.0 

12 

UNU4 

0 

LIN 

DO 

0 

100.0 

0.0 

If 

the  operator 

wants 

;  to 

enter  a  -25 

%  change 

f lowrate , 

he  wi 

1 1  enter  for 

parameter 

numbers 

2,9  return 

where  2  indicates  that  the  operator  wants  to  change  the 
VALUE  and  9  is  the  number  of  the  I/O-vector  (see  the  left 
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hand  margin  in  the  table  listing).  The  program  will  reply 
ENTER  NEW  VALUE,  0-10000: 
and  the  operator  would  enter 

2845  return  (i.e.  output  =  2845  mVolts)  which  will  set  the 
feed  flow  to  approximately  13.5  grams  per  second  (see 
explanation  in  next  section). 

3. 5. 6. 5  Feed  flow  rate  changes 

The  operator  can  set  the  set  point  for  the  feed  flow 
controller  by  manipulating  the  VALUE  in  the  I/O-vector 
label  led  FEED. 

Say  the  controller  has  a  bias  of  -0.25  grams/sec,  i.e.  if  we 
send  a  signal  corresponding  to  18  grams/sec  to  the 
controller,  the  actual  flow  will  be  regulated  to  17.75 
grams/sec.  This  so  called  "bias"  is  in  reality  the  sum  of 
biases  of  the  instruments  between  the  LSI- 11  and  the 
pneumatic  controller  (analog  output  module  in  LSI-11, 

Current  Output  Station,  I /P-conver ter  and  pneumatic  analog 
control ler ) . 

Given  this  -0.25  grams/sec  bias,  we  want  to  calculate  the 
output  voltage  for  the  feed  flow  so  that  the  actual  feed 
will  go  to  18  grams/sec. 

For  a  flow  rate  the  calculation  is: 

2 

Adesired  value  -  bias\ 

I  -  *  1 0000  mVol ts 

V  SPAN  ) 

Thus,  to  get  a  feed  flowrate  of  18  grams/sec,  the  output  has 


« 


71 


to  be 


18.0  +  0.25 

SPAN 


*  10000 


5013  mVolts 


3. 5. 6. 6  Event  scheduling 

To  enable  the  operator  to  schedule  events  ahead  of 
time,  which  facilitates  column  operation  without  the 
presence  of  an  operator,  an  event  scheduler  has  been 
incorporated  in  the  background  program.  This  basically 
consists  of  the  subroutine  SCHED  which  is  being  called  by 
subroutine  SERVE,  called  by  the  mainline  background  program. 
The  event  scheduling  requires  the  operator  to  enter 
parameter  changes  in  a  file  instead  of  on  the  system 
console.  The  filename  is  EVENT.DAT  in  floppy  drive  0. 
Theoretically,  an  operator  can  schedule  weeks  or  even  months 
ahead,  depending  on  the  size  of  the  event  file. 

Every  scheduled  event  consists  of  three  lines  of  code  in  the 
event  file: 

1.  The  first  line  contains  the  time  at  which  the  event  is 
to  take  place.  We  use  a  24  hour  clock  with  the  format 
hrs:min,  e.g.  1  am  is  01:00  and  4.15  pm  is  16:15.  This 
format  should  always  be  used  in  the  event  file. 

2.  The  second  line  contains  the  address  of  the  parameter 
that  is  to  be  changed.  The  format  for  this  line  is: 
table  number,  parameter  number  1,  parameter  number  2 


i 
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e.g.  the  feed  flowrate  setting  is  parameter  2,9  in  the 
I/O-table,  so  the  address  of  the  feed  flowrate  setting 
i  s 

2,2,9  table  number  2,  2nd  entry  in  the  I/O-vector, 

I/O -vector  number  9. 

The  set  points  are  the  third  parameter  in  the  CP-table, 
so  their  address  is 
3,3,0  table  number  3,  3rd  entry 
3.  The  third  line  contains  the  new  value  of  the  parameter 
that  is  to  be  changed,  e.g. 

2845  for  a  feed  flowrate 

95.0,4.0  for  the  set  points  at  95%  (top)  and  4%  (bottom) 

Example 

Say  it  is  4  pm  and  the  distillation  column  is  operating 
at  steady  state  (feed  flowrate  18  grams/sec,  top 
composition  95%,  bottom  composition  5%).  The  operator 
wants  to  schedule  the  following  events: 

-at  6  pm  the  feed  flowrate  should  drop  to  13.5  grams/sec 
-at  1  am  (next  day)  the  set  point  for  the  bottom 
composition  should  drop  to  4%  (wt%  methanol). 

EVENT . DAT 

18:00 

2,2,9, 

2845 


01:00 


,  ty  '  I  c?'"  "  *  ..  ^ 
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3,3,0, 

95.0,4.0, 

00:00 

NOTE:  After  the  last  scheduled  event  there  should  always  be 
a  line  00:00  to  signal  the  scheduler  that  it  should  stop 
reading  from  the  file.  Without  this  line  an  input  error 
wi 1 1  resul t . 

3.5.7  Modifying  the  Programs 

If  a  user  wants  to  make  modifications  in  the  programs 
to  accomodate  his  own  type  of  control,  he  will  have  to  write 
his  own  control  subroutine  (C0NTR).  All  other  changes  in  the 
subroutines  are  made  necessary  by  the  fact  that  a  user  may 
want  to  use  different  controller  parameters,  ergo  the  only 
changes  that  have  to  be  made  are  in  the  subroutines  that  use 
the  controller  parameters:  to  list  them  (LISCP),  to  change 
them  (CHANCP)  and  to  read  their  initial  values  from  a  file 
(mainline  foreground  and  background  programs). 

We  will  look  at  what  exactly  has  to  be  done: 

a.  A  user  will  probably  have  different  parameters  in 

his  control  algorithm.  These  parameters  will  have  to 
be  declared  in  the  subroutine  and  be  accomodated  in 
the  CP-table,  i.e.  the  parameters  have  to  be 
EQUIVALENCEd  to  words  in  the  CP-table. 

In  the  user  written  subroutine  C0NTR  the  only 
requirements  are: 

1)  the  COMMON  blocks  BLK 1  (I/0-table)  and  BLK2 
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(CP-table)  have  to  be  inserted; 

2)  if  the  input  values  for  the  flow  controllers  are 
calculated  in  engineering  units,  the  values  have 
to  be  converted  to  machine  units  (mVolts)  with  a 
CALL  INGE  (see  listing  of  CONTR  and  INGE  in 
appendix  A  ) 

The  user  can  put  his  parameters  in  the  CP-table  with 
an  EQUIVALENCE  statement,  similar  to  the  way  it  was 
done  in  the  standard  control  algorithm  (see  listing 
of  CONTR  in  appendix  A  ) 

b.  If  the  user  wants  to  have  the  same  flexibility  and 
convenience  in  listing  and  changing  his  own 
controller  parameters,  some  changes  will  have  to  be 
made  in  subroutines  LISCP  and  CHANCP  in  the 
subroutine  library  SUBBA.FOR. 

1)  LISCP 

The  only  requirement  is  that  COMMON  /BLK2/CPTAB 
be  declared.  The  rest  of  the  subroutine  consists 
of  an  EQUIVALENCE  statement  and  a  WRITE 
statement.  The  user  has  to  EQUIVALENCE  his  own 
parameters  to  the  CP-table  and  write  them  in  any 
format  he  1 i Kes . 

2)  CHANCP 

In  CHANCP  we  want  to  read  new  values  for 
controller  parameters.  The  user  again  has  to 
declare  COMMON  /BLK2/CPTAB  and  EQUIVALENCE  it  to 
his  own  parameters,  which  have  to  be  numbered 
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and  refered  to  by  the  numbers  chosen  by  the 
user  . 

The  simplest  way  of  modifying  CHANCP  is  to 
replace  the  EQUIVALENCE  and  READ  statements  in 
the  standard  subroutine  (see  listing  of  CHANCP 
in  appendix  A  ) 

c.  If  the  initialization  sections  of  the  foreground  and 
background  programs  read  all  initial  values  for  the 
parameters  from  a  file  INIT.DAT,  the  user  has  to 
make  the  following  changes  in  the  initialization 
sections  of  both  foreground  and  background  programs: 

1)  Replace  the  EQUIVALENCE  statement  for  the 
CP-table 

2)  Change  the  READ  statement  which  reads  data  from 
file  INIT.DAT  to  read  the  user's  parameters. 

3)  Of  course  the  user  also  has  to  replace  the 
present  parameter  values  in  file  INIT.DAT  by  the 
values  for  his  own  parameters. 

d.  If  the  user  wants  to  include  other  parameters  in  the 
data  acquisition,  he  will  have  to  include  those 
parameters  in  the  WRITE  statement  in  the  foreground 
program  which  writes  data  to  disk  file.  (See  listing 
of  F0R2.F0R  in  appendix  A) 


The  amount  of  effort  required  for  modifying  the  control 
subroutine  depends  on  the  complexity  of  the  user  s 
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algorithm,  which  replaces  the  current  robust  algorithm.  All 
other  changes  should  require  not  more  than  four  hours  of 
cod i ng . 

If  the  user  needs  more  controller  parameters  than  will 
fit  into  the  present  CP-table  (40  words),  the  user  will  have 
to  increase  the  size  of  the  CP-table.  The  complication  that 
arises  here  is  the  parameter  transmission  between  foreground 
and  background  programs.  The  user  will  have  to  increase  the 
size  of  the  transmission  buffers  as  well.  See  the  program 
listings  (appendix  A)  and  [24]. 


3.6  Evaluation  of  the  Control  System 

The  computer  system  we  used  to  implement  the  control 
algorithm  on  proved  to  be  easy  to  use  and  reliable.  The  fact 
that  we  have  a  dedicated  single  user  system  makes  it  in  some 
respects  more  flexible  and  dependable,  especially  where  the 
communication  with  the  GC  computer  is  concerned  since  the 
data  transmission  is  synchronous  (the  GC  does  not  wait  for 
an  LSI  signa 1 ) . 

After  two  weeks  of  testing  and  modifying,  all 
experimental  runs  on  the  distillation  column  were  done 
during  7  weeks  of  continuous  operation.  During  those  7  weeks 
we  experienced  9  system  failures,  3  caused  by  untraceable 
system  errors,  4  caused  by  GC  failures  and  2  caused  by 
erroneous  operator  entries.  The  system  failures  just  caused 
the  control  values  to  "freeze" ,  which  didn'  t  lead  to 
instability  or  a  crash  of  the  distillation  process. 


< 


. 


n  ’ 


4.  Experimental  Robust  Control  of  the  Distillation  Column 


4.1  Introduction 

Using  the  control  system  that  was  developed  in  the 
previous  part  of  this  thesis  we  can  now  evaluate  the 
performance  of  Davison's  robust  controller  on  a  real-life 
plant:  the  department's  pilot  plant  distillation  column. 

We  will  have  a  brief  look  at  distillation  column 
control  in  general  and  describe  the  set  up  for  two  quality 
energy  balance  control  which  we  used  for  the  implementation 
of  robust  control . 

Then  we  will  describe  the  distillation  column  that  was 
used  to  implement  the  control  on.  An  extensive  description 
of  the  control  hardware  can  be  found  in  part  2  of  this 
thesi s . 

Before  we  implement  a  controller,  we  have  to  perform 
the  steady  state  experiments  as  outlined  in  the  theory 
(section  2.2).  These  experiments  will  show  that  a  cascaded 
control  arrangement  with  proportional  control  in  the  slave 
loop  is  necessary. 

Then  the  design  and  implementation  of  the  robust 
controller  is  straightforward  and  we  can  evaluate  its 
performance  and  compare  it  to  conventional  multi  loop  PID 
control lers . 

One  of  the  problems  encountered  in  controlling  the 
distillation  column  was  the  long  time  delay  in  the  bottom 
loop.  The  conventional  strategy  for  dealing  with  time  delays 
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is  the  use  of  derivative  action  in  controllers.  We  tried  to 
improve  the  robust  controller  performance  by  including 
derivative  action  in  the  slave  loop  of  the  cascaded  robust 
controller.  We  will  see  that  robust  control  in  a  cascaded 
arrangement  indeed  does  offer  some  advantages  over 
conventional  P I ( D )  control. 


4.2  Distillation  Column  Control 

Separation  of  light  and  heavy  components  in  a  mixture 
is  one  of  the  most  important  operations  in  chemical 
industry,  especially  in  hydrocarbon  fuel  processing.  As  a 
direct  result  continuous  distillation  units  are  considered 
to  be  the  work  horses  of  most  chemical  plants  and  certainly 
of  oil  refineries  [27]. 

The  technique  of  distillation  is  well  developed  and 
major  gains  through  better  distillation  column  design  and 
construction  are  not  expected.  Nevertheless  the  rapidly 
increasing  energy  costs  have  provided  a  strong  incentive  to 
try  to  increase  the  efficiency  of  distillation  columns.  This 
has  led  to  an  increased  interest  in  control  [28*30]  and  the 
introduction  of  more  and  more  advanced  control  strategies, 
all  aiming  at  reducing  the  operating  costs. 

Especially  illustrative  is  optimal  control  where  the  control 
action  is  directly  aimed  at  minimizing  a  "cost  function". 
Most  control  stategies  however,  have  as  primary  objective 
operating  the  column  closer  to  the  specification  limits  with 
as  justification  an  expected  reduction  in  operating  costs. 
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Traditionally,  several  basic  control  strategies  have 
been  used,  depending  on  product,  objectives,  etc.  [27]: 

1.  Mass  balance  control.  Here  we  Keep  two  of  the  three  flow 
rates  (or  ratios)  constant.  This  method  can  be  used  when 
the  feed  composition  is  practically  constant  and  there 
is  one  product  which  is  much  more  important  than  the 
other . 

2.  Temper  at ure- and -pressure  control.  The  principle  is  to 
Keep  temperature  and  pressure  on  some  tray  in  the  column 
at  constant  values.  For  binary  mixtures  at  bubble  point 
the  composition  will  then  be  constant  too.  For  mixtures 
consisting  of  homologous  series  of  hydrocarbons  this 
method  generally  worKs  too,  provided  one  does  not  try  to 
fractionate  closely  boiling  components. 

3.  Single  quality  control.  We  measure  a  local  composition 
directly,  independent  of  the  pressure.  The  difference 
with  2)  lies  mainly  in  the  different  interations  between 
the  various  loops. 

All  the  above  methods  try  to  Keep  a  distillation  column 
close  to  an  ideal  operating  point.  If  we  would  want  to 
change  the  specifications  for,  say,  the  bottom  product,  but 
not  for  the  top  product,  the  change  in  operating  point  would 
not  be  straight  forward.  The  situation  is  more  clear  if  we 
use  two-quality-control,  with  composition  analyses  of  the 
top  and  bottom  products.  With  conventional  multi  loop  PID 
controllers,  we  have  two  direct  control  loops.  In  the  case 
of  the  departments  pilot  plant  column  (see  section  4.3) 
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these  loops  are: 

top  loop  -  in  which  the  top  product  composition  is  measured 
(in  the  condenser),  and 
the  reflux  flow  rate  is  manipulated, 
bottom  loop  -  the  bottom  product  composition  is  measured 
(sample  from  the  reboiler),  and 
the  reboiler  steam  flow  rate  is  manipulated. 

One  quality  controller  (bottom)  in  fact  adjusts  the  vapour 
flow  in  the  column,  the  other  (top)  the  liquid  flow. 

If  we  look  at  the  manipulated  variables  in  the  column, 
this  control  scheme  can  be  called  an  energy  balance  control 
scheme  because  we  manipulate  the  energy  flows  to  the  column 
(steam  and  reflux).  The  opposite  would  be  a  material  balance 
control  scheme  where  we  manipulate  product  flows  to  regulate 
compos i t ions . 

Tests  performed  by  Lieuson  [32]  showed  that  on  this 
column  for  the  operating  conditions  we  used,  an  energy 
balance  scheme  gives  better  results. 

Past  experiences  with  the  department's  distillation 
column  [30-40]  have  shown  that  controlling  this  column 
constitutes  a  realistic  test  for  any  control  algorithm, 
where  non  1 i near i t i es ,  time  delays  and  interactions  pose  the 
problems  that  are  to  be  overcome. 
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4.3  Distillation  Column  Equipment 

The  22.8  cm  diameter,  8  plate,  pilot  scale  distillation 
column  used  in  this  study  has  been  used  in  a  number  of 
previous  control  studies  [30-40].  Each  plate  at  a  spacing  of 
30.5  cm  is  fitted  with  4  bubble  caps.  Operating  conditions 
are  adjusted  to  provide  a  separation  of  the  50  weight 
percent  methanol -water  mixture  into  top  and  bottom 
compositions  of  about  95  and  5  weight  percent  methanol 
respectively.  A  detailed  summary  of  typical  steady  state 
operating  conditions  are  given  in  table  6. 

A  schema t i c  diagram  of  the  column,  which  is  completely 
instrumented  with  conventional  industrial  control  valves, 
controllers  and  transmitters  is  shown  in  figure  14.  As  can 
be  appreciated,  only  the  principal  control  instrumentation 
pertinent  to  the  objective  of  this  control  study  has  been 
shown.  Since  the  objective  is  to  maintain  top  and  bottom 
compositions  at  their  specified  values,  both  streams  are 
equipped  with  instrumentation  for  composition  analysis.  Top 
composition,  although  continuously  monitored  using  a 
temperature  compensated  in-line  capacitance  probe,  is 
sampled  at  the  same  rate  as  the  bottom  composition  signal. 
Bottom  composition  is  measured  using  a  HP  5720A  gas 
chromatograph  (GC).  Liquid  sampling  of  the  bottom  product 
stream  is  under  the  control  of  one  of  three  HP  1000  s  of  the 
Department's  distributed  computer  network  of  mini  and 
micro-computers.  Analysis  of  the  chromatogram,  on  a  180 
seconds  cycle,  is  also  performed  by  this  node  of  the 
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Table  6 . 

Typical  steady  state  operating  conditions 


Feed  flow  rate . 18.0  g/s 

Reflux  flow  rate .  8.8  g/s 

Steam  flow  rate . 13.6  g/s 

Top  product  flow  rate....  9.0  g/s 

Bottom  product  flow  rate  9.0  g/s 

Feed  composition . 50.0  %  (weight  %  methanol) 

Top  composition . 95.0  %  (weight  %  methanol) 

Bottom  composition . 5.0  %  (weight  %  methanol) 


distributed  network  wich  then  transmits  a  signal  to  a 
' local7  node  of  the  network,  a  DEC  11/03  microcomputer.  The 
control  sample  interval  for  both  loops  is  equal  to  3 
minutes,  as  a  unit  integer  multiple  of  the  3  minute 
measurement  delay  time  that  occurs  in  the  GC  analysis.  The 
schematic  configuration  of  the  pilot  plant  units  and  the 
distributed  computer  network  is  shown  in  figure  9  (page  48). 
As  can  be  seen  from  the  diagram,  implementation  of  the 
control  algorithm  and  data  acquisition,  except  for  the  GC 
composition  measurement  ,  is  performed  using  the  DEC  11/03. 
Operator  communication  to  the  DEC  11/03  is  by  means  of  a  CRT 
terminal.  See  section  3.4  for  a  more  detailed  description  of 
the  system. 


4.4  Experiments  for  Controller  Design 

As  was  outlined  in  the  description  of  the  theory,  the 
robust  controller  design  technique  was  intended  for  use  on  a 
linear  system.  In  section  2.3  we  have  shown  that  a  nonlinear 
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Figure  14.  Pilot  plant  distillation  column  with 
instrumentation 
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plant  may  yield  different  controllers  for  the  same  operating 
points,  depending  on  whether  we  use  positive  or  negative 
inputs.  This  is  also  the  case  for  the  distillation  column 
under  consideration. 

For  positive  step  inputs  (+10%  reflux  flow  rate  and  +8% 
steam  flow  rate)  we  find: 


'0.826 

-1 . 400~ 

u;  ]■'  = 

'-7.785 

2.477' 

(37) 

.3.000 

-4.400, 

.-5.308 

1 .462. 

For  negative  step  inputs  (-10%  reflux  flow  rate  and  -8% 
steam  flow  rate)  we  find: 


1 

'1.077 

-0.539“ 

i  —  / 

'1.149 

o 

o 

1 

I T  ,  1  = 

n:  ]  = 

2.366 

-6. 1 54  _ 

0.442 

-0.201, 

Appendix  B  plots  19  and  20  show  responses  for  positive  and 
negative  changes  in  reflux  flow  rate. 

We  see  that  (37)  and  (38)  have  opposite  signs,  so  we 
cannot  construct  a  controller  for  this  operating  point.  The 
differences  in  (37)  and  (38)  are  such  that  averaging  does 
not  make  sense  since  we  cannot  even  determine  the  signs  in 
the  controller  matrix  from  above  results.  Therefore  we  first 
apply  multi  loop  proportional  feedback  control  to  the 
distillation  column,  as  was  suggested  in  section  2.4  . 

The  primary  purpose  of  the  P-control  is  to  "linearize" 
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the  plant,  i.e.  to  give  as  much  as  possible  the  same 
responses  for  positive  and  negative  set  point  changes.  We 
want  to  avoid  the  behaviour  as  outlined  in  1.3  where  we 
showed  that  the  robust  design  technique  on  a  nonlinear  plant 
can  give  us  two  totally  different  controllers  for  the  same 
set  of  operating  conditions. 

The  choice  of  Kp  is  not  very  important  and  only  a  few 
test  runs  are  necessary  to  obtain  acceptable  values.  Since 
values  for  PID-control lers  were  available,  we  started  with 
the  gains  from  the  PID-control lers  (3.571  top  and  -0.576 
bottom).  Plot  21  (appendix  B)  shows  the  responses  to  +20% 
change  in  feed  flow  rate.  The  response  is  fast  with  a  little 
bit  of  overshoot  in  top  and  bottom  composition.  A  full 
series  of  experiments  as  outlined  in  [1]  (see  section  2.2.2) 
for  both  positive  and  negative  inputs  and  disturbances  shows 
that  the  objective  of  "plant  linearization"  is  fairly  well 
sat i sf i ed : 

positive  changes: 


1.506  -0.032 


r-0.0549 


IT  I 


D 


(39) 


-1.649  1.359 


0.7985 


negative  changes: 


1.306  -0.069 


-0.0216 


D 


0.5696 


(40) 


-0.979  1.488 
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average : 


1  -1 

'  1.406 

i 

LO 

O 

O 

l 

i 

* -0 . 038  ' 

IT  1  = 

Di  = 

1 

.-1.314 

1 .424, 

.  0.684, 

Because  the  response  shows  some  overshoot  which  might 
lead  to  oscillatory  behaviour  in  a  robust  controlled 
situation,  we  repeat  the  experiments  with  lower  gains.  We 
use  70%  of  the  previous  gains  (2.5  top  loop,  -0.4  bottom 
loop).  Because  the  offset  will  be  larger  we  decrease  the 
step  size  of  the  inputs.  Plots  22  and  23  (appendix  B)  show 
the  responses  to  +10%  and  -10%  changes  in  the  feed  flow 
rates.  This  time  we  have  fast  and  smooth  responses. 
Calculation  of  [ T J  ]  -matrices  after  conducting  all 
experiments  shows  that  the  linearizing  effect  of  P-control 
is  still  satisfactory: 


positive  changes: 


1  .370 

-0.077' 

i 

"-0 . 026  " 

= 

D,  = 

.-1.841 

1  .423. 

i 

.  0.935. 

(42) 


negative  changes: 


1.581  0.111 
-2.096  1.988 


■'-0.042" 

0.744 


(43) 
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average : 


1  .476 

0.017' 

■ 

0 . 034  ' 

- 

D  = 

.-1  .969 

1  .706. 

1 

.  0.840  . 

(44) 


To  round  off  the  experiments  we  repeat  the  same 
procedure  for  gain  values  of  50%  of  the  first  case  (1.785 
top,  -0.288  bottom).  Plot  24  (appendix  B)  shows  the 
responses  to  +10%  change  in  feed  flow  rate.  The  transients 
are  smooth  and  slightly  slower  than  in  the  previous  case. 
There  is  still  a  reasonable  amount  of  linearization: 


positive  changes: 


i  —  | 

'2.770 

-0.469" 

i 

"-0 . 054 " 

it;i  1  = 

.-5.994 

2.983. 

D,= 

.  0.831 . 
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negative  changes: 


1.974  -0.116' 

-3.212  1.922 


“-0.018 
.  1.057 
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average : 
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We  found  that  the  second  case  with  gains  2.5  for  top 
loop  and  -0.4  for  bottom  loop  gave  an  acceptable  amount  of 
plant  linearization  combined  with  smooth  and  fast  responses 
that  will  probably  provide  the  best  basis  for  the  next  step, 
the  implementation  and  tuning  of  a  robust  controller. 

4.5  Robust  Feedback  Feedforward  Control 

We  are  most  interested  in  regulatory  control,  i.e.  the 
ability  of  the  controller  to  keep  the  output  values  (top  and 
bottom  product  compositions)  as  close  as  possible  to  the 
desired  values  (set  points)  in  the  face  of  disturbances  and 
parameter  variations. 

Traditionally,  controllers  for  the  distillation  column 
have  been  tuned  for  feed  flow  disturbances  while  the  set 
points  were  being  kept  constant.  The  robust  controllers  that 
will  be  described  in  this  and  the  following  sections  will  be 
tuned  for  a  -25%  feed  flow  rate  disturbance.  The  PID 
controllers  that  were  available  for  comparison  were  also 
tuned  for  this  disturbance. 

The  controller  we  used  was  (see  section  4.4): 


o 


where : 


f- 105  -9.80 


0.00377 
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The  controller  was  tuned  for  a  -25%  feed  flow  disturbance, 
both  with  and  without  feed  forward  action  (i.e.  D,'  =  0). 

The  tuning  was  done  by  manipulating  epsilon  and  calculating 
the  Integral  Absolute  Errors  (IAE)  for  top  and  bottom 
product  composition  from  the  responses.  See  table  7  for  the 
results.  (Calculation  over  1.5  hrs.  =  30  samples  was  carried 
out,  which  was  sufficient  for  the  transient  to  die  out) 

I f  we  look  at  the  sum  of  IAEs  for  top  and  bottom 
product  compositions,  we  find  that  the  best  controller 
setting  for  robust  control  without  feedforward  is  epsilon  = 
0.15.  For  robust  feedback  feedforward  control  we  find 
epsilon  =  0.10  the  best  setting. 

A  full  series  of  runs,  including  feed  flow  disturbances 
and  set  point  changes  for  the  best  controller  settings  gave 
results  as  outlined  in  table  8. 

Appendix  B  plots  25  to  30  show  responses  for  feed  flow  and 
set  point  changes  under  robust  control. 

For  robust  control  without  feedforward  we  find  the  IAE 
values  outlined  in  table  9. 


4.6  Comparison  of  Robust  and  PID  controllers 

In  this  section  we  will  look  at  the  robust  controller's 
performance  relative  to  conventional  multi  loop  PID 
control lers . 

A  PID  controller  uses  a  control  algorithm: 


u  =  Kp*E  +  Kj  *  2  E  +  K0*(yo-y) 


(49) 
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Table  7. 

Tuning  of  robust  feedback  feedforward  controller 
IAE  IAE 

wi thout  feedforward  wi th  feedforward 


epsi Ion 

top 

bottom 

top 

bottom 

0.05 

6.8 

58.8 

5.7 

4.9 

0.10 

4.4 

35.4 

3.0 

6.2 

0.15 

3.5 

30.7 

2.4 

9.0 

0.20 

3.9 

34.4 

2.2 

10.5 

0.25 

4.9 

40. 1 

2.2 

11.0 

Table  8. 

Load  disturbances  and  set  point  changes  for 

robust  feedback  feedforward  control 

IAE  IAE 

return  to 


steady  state 


CHANGE 

top 

bottom 

top 

bot  tom 

+25%  feed  flow  rate 

4.5 

9.3 

4.3 

4.2 

-25%  feed  flow  rate 

3.0 

6.2 

3.4 

5.  1 

+1%  top  composition 

3.  1 

5.7 

3.4 

7.2 

-1%  top  composition 

2.2 

4.3 

1.9 

4.3 

+1%  bottom  composition 

0.8 

10.2 

0.9 

7.4 

-1%  bottom  composition 

0.8 

7.3 

1  .0 

8.0 

. 

' 
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Table  9. 


Load  disturbances  and  set  point  changes  for 

robust  feedback  control 


CHANGE 


+25%  feed  flow  rate 
-25%  feed  flow  rate 

+1%  top  composition 
-1%  top  composition 
+1%  bottom  composition 
-1%  bottom  composition 


IAE  IAE 

return  to 


steady  state 


top 

bottom 

top 

bottom 

4.3 

29.8 

3.7 

28.3 

4.0 

28.2 

3.8 

24.8 

2.  1 

5.8 

2.7 

6.  1 

2.9 

6.0 

2.0 

4.7 

1  .3 

8.9 

0.9 

7.4 

0.8 

7.3 

1  .0 

7.8 

where  u  - 


K 


o 


control  value 

proportional  gains 

integral  gains  =  (KP*ISEC)/TI 

where  ISEC  =  sample  interval  in  seconds 
Tx  =  integral  time  in  seconds 
derivative  gains  =  (KP*T^)/ISEC 

where  T^  =  derivative  time  in  seconds 


E  =  error  =  yvc|  -  y 

y0  =  previous  value  of  the  output  y 


The  controller  we  used  for  comparison  was  tuned  for  a  -25% 
feed  flow  disturbance.  The  controller  settings  are: 

Top  loop  :  y  =  top  product  composition 
u  =  reflux  flow  rate 
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Kp=  3.571,  Tj.  =  695.0  secs.,  Tj  =  41.6  secs. 

Bottom  loop  :  y  =  bottom  product  composition 

u  =  reboiler  steam  flow  rate 
Kp  =  -0.576,  Tj  =  771.0  secs.,  =  39.5  secs. 

The  performance  comparison  for  various  disturbances  and  set 
point  changes  is  outlined  in  table  10  .  The  PI  controller  is 
the  same  as  PID  but  with  the  derivative  action  set  to  zero. 
Appendix  B  plots  31  and  32  show  -25%  feed  flow  rate 
responses  for  PI  and  PID  control. 

The  above  results  show  that  PI  and  PID  control  is  slightly 
better  than  robust  control.  The  robust  controller  however 
was  designed  with  feedforward  action  in  mind.  It  is  obvious 
that  the  performance  of  the  robust  feedback  feedforward 
controller  for  load  disturbances  will  be  better  than  a 
normal  PID  controller,  due  to  the  feedforward  action.  To 
make  the  comparison  more  realistic,  we  added  feedforward 
action  to  the  PID  controller  so  that  the  PID  +  feedforward 
control  algorithm  is: 

u  =  Kp  *E  +  Kx  *  5  E  +  K  0  *  ( y„  -  y )  +  Kp  *  [  T,'  ]  '*Dt'  *FF  (50) 

FF=  Feed  flow  rate 

See  equation  (41)  for  [  T  J  ]  '  and  D,'  . 

Note  that  the  [ T | ]"'  and  D ‘  matrices  were  designed  for  the 
values  of  Kp  used  in  the  PID  controller  (see  section  4.4) 
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Table  10. 

Load  disturbances  and  set  point  changes  for 

robust  feedback  control  vs.  PI  and  PID  control 

IAE 


Robust 

feedback 

PI 

PID 

CHANGE 

top 

bot  tom 

top 

bot  tom 

top 

bot  tom 

+25%  feed  flow 

4.3 

29.8 

2.6 

25.9 

2.8 

27.3 

-25%  feed  flow 

4.0 

28.2 

2.3 

23.6 

2.  1 

21  .7 

+1%  top  comp. 

2.  1 

5.8 

3.0 

11.6 

3.  1 

10.0 

-1%  top  comp. 

2.9 

6.0 

2.4 

6.9 

3.0 

8.4 

+1%  bottom  comp. 

1  .3 

8.9 

0.8 

6.  1 

0.7 

4.4 

-1%  bottom  comp. 

0.8 

7.3 

0.8 

5.2 

1  .  1 

5 . 6 

15.4  86.0 

11.9  79.3 

12.8 

77.4 

101.4 

91  .2 

90 

.2 

The  performance  comparison  is  given  in  table  11.  Appendix  B 
plots  33  to  35  show  responses  to  feed  flow  and  set  point 
changes  for  P ID-p 1  us- feedforward  control. 

We  see  that  the  performance  of  the  robust  controller  is 
roughly  the  same  as  the  PI  and  PID  performance.  For  set 
point  changes  we  find  that  the  robust  controller  is  better 
for  changes  in  top  composition,  worse  for  bottom  composition 
changes.  This  can  be  explained  from  the  fact  that  there  is  a 
strong  interaction  between  the  steam  flow  rate  and  the  top 
product  composition.  The  multivariable  control ler  accounts 
for  this  interaction  and  thereby  minimizes  the  undesirable 
change  in  bottom  composition.  In  case  of  a  set  point  change 
in  bottom  composition,  the  robust  control ler  basically 
behaves  as  a  conventional  multi “loop  controller  since  there 
is  hardly  a  ripple  in  the  top  composition.  Here  the  robust 
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Table  11. 

Load  disturbances  and  set  point  changes  for 

robust  feedback  feedforward  control  vs.  PI  and 

PID  control 


Robust  feedback 

feedforward 

IAE 

PI 

+f eedforward 

PID 

+feedforward 

CHANGE 

top 

bot  tom 

top 

bottom 

top 

bottom 

+25%  feed  flow 

4.5 

9.3 

1.8 

9.5 

1.5 

9.  1 

-25%  feed  flow 

3.0 

6.2 

4.3 

23.7 

1  .7 

12.6 

+ 1  %  top  comp . 

3.  1 

5.7 

3.0 

11.6 

3.  1 

10.0 

-  1%  top  comp . 

2.2 

4.3 

2.4 

6.9 

3.0 

8.4 

+1%  bottom  comp. 

0.8 

10.2 

0.8 

6.1 

0.7 

4.4 

-1%  bottom  comp. 

0.8 

7.3 

0.8 

5.2 

1  .  1 

5.6 

14.4 

43.0 

10.5 

62.4 

11.1  50.1 

57 

.4 

72, 

.9 

61.2 

controller  loses  its  multivariable  advantage  over  the  PID 
control ler . 

If  we  look  at  the  sum  of  IAEs,  we  see  that  PI  and  PID 
controllers  do  a  slightly  better  job  controlling  the  top 
loop,  but  the  robust  controller  gives  better  control  for  the 
bottom  loop.  This  means  that  the  robust  controller's 
advantage  in  i ncorporat i ng  interaction  in  the  control 
algorithm  outweighs  the  tuning  advantages  and  derivative 
action  in  the  PID  controllers.  The  bottom  loop  is  considered 
to  be  more  difficult  to  control  than  the  top  loop  because  of 
the  long  time  delays.  The  robust  control ler  has  a  slight 
advantage  over  PID  controllers  since  the  IAE  presents  a 
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measure  of  the  amount  of  off-spec  product  produced  by  the 
distillation  column.  This  is  less  than  for  PID  control. 


4.7  Robust  Control  plus  Derivative  Action 

The  major  control  problem  in  our  system  is  posed  by  the 
long  time  delay  in  the  bottom  loop  (7  minutes).  Because  of 
its  phase  lead  properties  derivative  action  has  been 
proposed  as  medicine  for  time  delays  (phase  lag)  in 
conventional  control  schemes.  It  may  be  expected  that  adding 
derivative  action  to  the  slave  loop  of  the  robust  controller 
will  give  us  some  improvement  in  the  control.  Since 
derivative  action  operates  on  the  system  outputs  only,  it  is 
easy  to  add  to  the  control  algorithm.  Also,  derivative 
action  does  not  affect  the  steady  state  behaviour  of  a 
controller,  therefore  we  do  not  have  to  redesign  the  robust 
controller  if  we  want  to  use  PD  control  in  the  slave  loop 
instead  of  P  control. 

Since  we  want  to  Keep  as  much  as  possible  the  tuning 
advantage  of  the  robust  controller,  the  derivative  action 
was  implemented  as  "add-on"  feature  and  also  tuned  as  such, 
i.e.  the  best  robust  controller  tuning  from  the  previous 
data  was  taken  and  then  only  the  D  action  was  varied  in  the 
tuning  runs.  Again  all  tuning  runs  were  -25%  feed  flow 
disturbances.  The  tuning  results  are  listed  in  table  12 
As  we  see  the  effects  of  derivative  action  are  minimal .  The 
differences  in  IAEs  cannot  be  call ed  signi ficant .  Our  bes t 
choice  of  derivative  time  40  seconds  for  both  loops  is 
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Table  12. 

D-tuning  for  -25%  feed  flow  disturbance 

PD-robust  feedback  feedforward  control 


Derivative  time  IAE 


top  loop 

bottom  loop 

top 

bot  tom 

10 

10 

4.4 

6.3 

20 

20 

4.3 

7.0 

30 

30 

4.2 

5.3 

40 

40 

4.2 

5.0 

50 

50 

4.2 

6.  1 

10 

20 

3.8 

7.4 

20 

40 

4.2 

7.5 

40 

80 

4.2 

8.7 

80 

160 

3.7 

8.4 

160 

320 

4.7 

10.2 

somewhat  based  on  intuition  and  supported  by  the  fact  that 
this  is  the  same  amount  of  derivative  action  as  in  the  well 
tuned  PID  controller. 

We  also  did  some  derivative  tuning  for  PD-robust 
control  without  feedforward  action  but  the  differences 
between  the  runs  were  very  insignificant.  The  comparison  of 
the  overall  controller  performances  is  given  in  tables  13 
and  14.  Appendix  B  plots  36  to  41  show  the  column  responses 
to  feed  flow  disturbances  and  set  point  changes  for 
PD-robust  control . 

As  we  can  see,  the  overall  results  for  PD-Robust 
control  are  slightly  better  than  P-Robust  or  PID  control. 
The  control  improvement  lies  largely  in  the  bottom  loop 
(compare  the  total  IAE  figures  for  the  bottom  loops).  This 
was  to  be  expected  since  the  derivative  action  can  cope 
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Table  13. 

Load  disturbances  and  set  point  changes  for 

P-  and  PD-robust  feedback  control  vs.  PID 
control 


P- Robust 


CHANGE 

top 

bot  tom 

+25%  feed  flow 

4.3 

29.8 

-25%  feed  flow 

4.0 

28.2 

+1%  top  comp. 

2.  1 

5.8 

-  1%  top  comp . 

2.9 

6.0 

+1%  bottom  comp. 

1  .3 

8.9 

- 1%  bottom  comp. 

0.8 

7.3 

15.4 

86.0 

101.4 


IAE 

PD- 

Robust 

PID 

top 

bottom 

top 

bot  tom 

4.2 

28.0 

2.8 

27.3 

4.6 

33.4 

2.  1 

21  .7 

2.  1 

4.7 

3.  1 

10.0 

3.7 

3.9 

3.0 

8.4 

1  .3 

8.4 

0.7 

4.4 

1  .4 

7.2 

1  .  1 

5.6 

17.3 

85.6 

12.8 

77.4 

102. 

9 

90 

.2 

Table  14. 


Load  disturbances  and  set  point  changes  for 

P-  and  PD- Robus t+feedforward  control  vs. 
P I D+Feedforward  control 


P-Robust 

+f eedforward 

IAE 

PD-Robust 

+feedforward 

PID 

+f eedforward 

CHANGE 

top 

bot  tom 

top 

bot  tom 

top 

bottom 

+25%  feed  flow 

4.5 

9.3 

4.0 

2.7 

1.5 

9.  1 

-25%  feed  flow 

3.0 

6.2 

4.2 

5.0 

1  .7 

12.6 

+1%  top  comp. 

3.  1 

5.7 

2.5 

3.6 

3.  1 

10.0 

-1%  top  comp. 

2.2 

4.3 

3.8 

4.6 

3.0 

8.4 

+1%  bottom  comp. 

0.8 

10.2 

1  .7 

7.0 

0.7 

4.4 

-  1  %  bot  tom  comp . 

0.8 

7.3 

1  .3 

7.8 

1  .  1 

5 . 6 

14.4 

43.0 

17.5 

30.7 

11.1  50.1 

57 

.4 

48. 

2 

61  .2 

’ 
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better  with  the  long  time  delay  we  find  in  the  bottom  loop. 


' 


' 

'  | 


5 .  Cone  1 us i ons 


5.1  Controller  Evaluation 

The  evaluation  of  the  robust  controller  performance  can 
be  based  on  both  the  simulations  from  chapter  2  of  this 
thesis  and  the  experiments  in  chapter  4.  Throughout  this 
thesis  we  have  been  using  Integral  Absolute  Errors  (IAE)  for 
comparison  with  PID  controllers.  Before  we  start  drawing 
conclusions  for  the  "quality  of  the  control",  some  comments 
should  be  made  about  the  i nterpretat ion  of  IAEs. 

Firstly,  the  objective  of  control  is  to  Keep  the  output 
values  as  close  as  possible  to  the  desired  values  under  all 
circumstances.  For  a  constant  throughput,  the  IAE  represents 
a  measure  of  the  amount  of  off -spec  product  produced.  Since 
we  are  dealing  with  compositions,  we  can  mix  liquid  that 
exceeds  the  specification  limits  with  liquid  that  is  below 
the  specification  limit  to  obtain  an  end  product  that  is 
exactly  on-spec.  Similarly,  if  the  composition  of  one  of  the 
distillation  column  products  would  follow  a  sinusoidal 
oscillation  around  the  desired  value,  we  would  (on  the 
average)  produce  no  off-spec  product,  but  the  IAE  would  Keep 
on  increasing.  Therefore  it  is  sometimes  argued  that  in  the 
case  of  compositions  as  outputs,  an  Integral  Error  criteria 
would  be  better  than  an  Integral  Absolute  Error  criteria. 
Counter  arguments  are  that  oscillations  in  the  outputs  do 
represent  problems  for  downstream  units.  This  should  be 
reflected  in  the  performance  criteria.  Also  the  basic 
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objective  of  control  is  minimizing  the  operating  costs. 
Oscillations  in  product  quality  and  subsequent  blending  may 
produce  on-spec  product,  but  the  operation  represents  an 
entropy  increase  and  is  therefore  wasteful  of  energy  [29]. 
The  IAE  does  "penalize"  oscillations. 

Our  decision  to  use  IAE  was  based  on  these  considerations 
and  on  the  fact  that  IAE  results  for  PID  controllers  were 
available  for  comparison. 

Secondly,  in  the  case  of  simulations  the  IAE  values  are 
determined  by  system  dynamics  and  controller  performance 
only.  In  the  experimental  runs  we  also  have  to  contend  with 
noise,  parameter  drift  and  unmeasurable  disturbances. 
Therefore,  small  differences  in  IAE  values  between  PID  and 
Robust  controller  for  simulation  runs  are  the  direct  result 
of  controller  performance,  the  system  dynamics  being  the 
same.  But  in  the  experimental  runs,  minor  variations  in  IAE 
are  not  meaningful  1  for  the  controllers  because  they  might 
be  caused  by  external  disturbances,  noise,  etc. 

The  first  conclusion  that  should  be  drawn  is  that 
Davison's  robust  controller  as  such  cannot  deal  with  a 
system  exhibiting  the  non  1 i near i t i es  as  mentioned  in  this 
thesis.  Only  after  modifications  (the  addition  of  a  slave 
control  loop)  we  can  expect  reasonable  control ler 
performance . 

Looking  at  the  simulation  results,  we  see  that  robust 
feedback  control  performs  slightly  better  than  PI  control. 

If  we  add  feedforward  action  to  the  controllers,  the  gap 
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widens,  and  especially  for  the  more  important  load 
disturbances  we  see  a  significant  improvement.  Adding  up  all 
IAE  values  for  feed  flow  and  set  point  changes,  we  find  that 
the  total  IAE  for  robust  feedback  feedforward  control  is  19% 
less  than  for  PI  control.  Looking  at  feed  flow  disturbances 
only,  we  find  an  improvement  of  40%. 

If  we  add  derivative  control  action  to  both  Robust  and 
PI  controllers,  we  find  an  improvement  in  control.  Now  the 
total  IAE  improvement  is  29%  but  for  feed  disturbances  only 
an  improvement  of  3%,  ergo  the  improvement  lies  largely  in 
set  point  changes. 

For  the  experimental  runs  we  find  that,  compared  to 
PID-plus-feedforward  control,  P-robust  feedback  feedforward 
control  gives  a  6%  overall  improvement.  Including  derivative 
action  gives  control  21%  better  than  PID-plus-feedforward . 
Looking  at  feed  flow  disturbances  only,  we  find  P-robust 
control  8%  better  and  PD-robust  36%  better  than  PID  control. 

The  difference  in  improvements  between  simulations  and 
experimental  runs  is  caused  by  noise  and  unmeasurable 
disturbances  in  the  experimental  case  which  make  integral 
control  action  necessary.  In  the  simulations,  integral 
action  in  the  controller  actually  caused  a  deterioration  of 
control  performance  (see  section  2.5.4).  As  has  been  shown 
in  the  previous  chapters,  robust  multivariable  control  can 
offer  significant  advantages  over  conventional  control.  This 
conclusion  is  also  supported  by  [41],  although  this  paper 
does  not  indicate  any  measures  to  counter  nonlinearity 
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problems.  If  we  compare  all  controller  performances  that 
were  outlined  in  the  previous  sections,  we  find  that  the 
best  distillation  column  control  is  achieved  with  PD-Robust 
Feedback  Feedforward  control.  Although  the  improvements  over 
conventional  PID  control  do  not  seem  to  be  very  significant, 
one  should  keep  in  mind  that  distillation  column  control  is 
mostly  regulatory  control.  If  we  look  at  the  controller 
performances  for  feed  flow  disturbances  only,  we  see  that 
robust  controllers  in  our  case  do  offer  some  significant 
improvements  due  to  the  much  better  control  of  the  bottom 
loop . 

Also  one  should  consider  that  the  performance  ot  the 
PID  controller  is  upgraded  by  the  addition  of  feedforward 
control.  The  feedforward  controller  was  designed  with 
Davison's  technique,  so  part  of  the  credit  for  the  good 
showing  of  the  PID  controller  should  go  to  the  robust 
controller  design  technique. 

Although  it  does  not  show  up  in  plots  and  tables,  we 
should  also  mention  that  it  took  only  a  minimal  tuning 
effort  to  achieve  a  robust  controller  which  equals  or 
surpasses  the  performance  of  the  conventional  PID  scheme. 
Proper  tuning  of  a  PID  controller  is  a  very  time  consuming 
and  complicated  task  which  requires  a  high  degree  of 
expertise  from  the  control  engineer,  and  the  tuning  process 
would  probably  yield  much  more  off~spec  product. 
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5.2  Recommendations  for  future  work 

A  possible  improvement  can  be  made  in  the  control 
software.  Currently  the  control  programs  take  up  almost  all 
available  memory  in  the  LSI  so  that  the  system  would  not 
accomodate  a  much  larger  control  algorithm.  Some  subroutines 
can  be  made  more  space  efficient  without  affecting  the 
system  performance.  Especially  the  data  transmission  between 
foreground  and  background  programs  which  presently  require 
large  buffers  can  be  made  more  efficient  if  necessary. 

One  of  the  major  problems  in  the  control  is  the  long 
measurement  time  delay  in  the  bottom  loop  (3  minutes) 
because  of  the  GC  analysis.  Improvement  of  controller 
performance  could  be  achieved  by  reducing  this  time  delay. 
This  can  be  done  by  reducing  the  GC  analysis  from  2.5 
minutes  to  1.5  to  2  minutes.  Next,  one  could  choose  a  higher 
baud  rate  for  data  transmission  from  the  GC  computer  to  LSI 
which  can  reduce  the  report  time  from  approximately  30 
seconds  to  less  than  1  second.  Previously  this  was 
impossible  because  of  limitations  of  the  time-sharing 
control  computer  (HP  1000)  which  demanded  a  slow  data 
transmi ssion . 

An  even  bigger  improvent  could  be  achieved  by  using  the 
temperature  on  tray  1  as  an  indicator  for  bottom 
composition.  Since  we  have  a  binary  mixture  at  bubble  point, 
Gibb's  phase  rule  says  that  we  have  two  degrees  of  freedom. 
Since  this  is  an  atmospheric  distillation  column,  the 
pressure  is  fixed  and  by  measuring  the  temperature  we  can 
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determine  the  composition.  The  GC  analysis  could  still  be 
used  in  a  cascaded  arrangement  to  improve  accuracy. 

By  thus  enabling  a  greatly  increased  sample  rate  at  no 
measurement  delay,  a  major  control  improvement  may  be 
expected . 
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NAME:  F0R2 
DATE:  791114 


LOCATION:  FOR2.FOR 
AUTHOR:  G.W.M.  COPPUS 


ABSTRACT:  Mainline  foreground  program,  see  description  below 
CALLING  SEQUENCE: 


CALLED 

CALLS: 

BY: 

ASSIGN. . . . 

. . . . System 

Subrout i ne 

CLOSE . 

. . . . System 

Subrout i ne 

IQSET . 

Subrout i ne 

GTIM . 

Subrout i ne 

MTATCH . 

. . . . System 

Subrout i ne 

MTGET . 

. . . . System 

Subrout i ne 

MTSET . 

. . . . System 

Subrout i ne 

IRCVDF . 

, . . . System 

Subrout i ne 

CVTTIM . 

, . . . System 

Subrout i ne 

ISDATF . 

, . . . System 

Subrout i ne 

ISLEEP . 

,  .  .  . System 

Subrout i ne 

GLUE . 

...Assembler  Subroutine 

DOUT . 

. . . Assemb 1 

er  Subroutine 

GETBC . 

. . . For  tr an 

Subrout i ne 

TIMER . 

. . . Fortran 

Subrout i ne 

DATAC . 

. . . For  tran 

Subrout i ne 

CONTR . 

. . . Fortran 

Subrout i ne 

PUT . 

. . . For  tran 

Subrout i ne 

COMP . 

. . . Fortran 

Subrout i ne 

COMD . 

. . . For  tran 

Subrout i ne 
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DESCRIPTION: 

In ial i zat i on  sect i on 

After  the  declaration  statements,  file  INIT.DAT  on  drive 
DYO:  will  be  opened.  It  contains  all  initial  values  for 
variables  in  the  three  tables  (I/O-table,  CP-table, 

PP- table) . 

The  values  are  read  and  the  file  is  closed  again.  A 
separate  subroutine  (GLUE)  puts  the  real  values  of  SPAN 
and  ZERO  into  the  integer  arrays  of  the  I/O-table. 

IQSET  sets  aside  extra  space  for  queue  elements  needed 
for  timer  and  10  calls. 

The  GC  link  (see  description  of  subroutine  GETBC)  is  set 
up  as  a  terminal  link.  CALL  MTATCH  attaches  a  new 
terminal  to  the  system,  reserves  buffer  space  for  10, 
etc.  MTGET  returns  the  terminal  character i st ics  and 
MTSET  sets  the  new  characteristics  (e.g.  inhibit  echos, 
see  GETBC ) . 

The  GC  remote  start  switch  is  closed  by  default  at 
system  power  up.  The  program  sets  the  switch  to  its 
normal  (open)  position. 

The  program  opens  a  file  for  data  acquisition.  The 
filename  has  to  be  entered  by  the  operator. 

Main 1 ine  sect  ion 

The  program  checks  a  flag  to  see  if  there  were  any 
parameters  sent  by  the  background  program.  If  yes,  the 
new  parameters  are  put  in  the  appropriate  tables,  and  a 
new  "receive  parameter  request"  will  be  issued. 
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An  ITYPE  greater  or  equal  than  20  is  the  signal  that  the 
operator  wants  to  close  the  data  acquisition  file  and 
open  a  new  one . 

The  GC  report,  sent  by  the  HP  1000  GC  computer,  is  read 
in  subroutine  GETBC 

The  timer  waits  for  the  next  sample  time. 

The  GC  is  triggered  (by  a  digital  output  signal)  to  get 
a  new  sample. 

All  other  process  variables  are  read  by  subroutine 
DATAC. 

Control  settings  are  calculated  in  subroutine  CONTR. 

The  signals  to  the  valve  controllers  are  sent  by 
subroutine  PUT. 

Every  sample  time,  the  most  important  data  are  being 
sent  to  disk  file. 

If  the  background  program  is  running,  data  will  be  sent 
for  operator  presentation. 


, 
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NAME:  DOUT  LOCATION:  DOUT.MAC 

DATE:  800114  AUTHOR:  R.L.  BARTON 


ABSTRACT:  Digital  output  driver,  to  set  the  digital  outputs 
on  the  LSI- 11  to  the  desired  values. 

CALLING  SEQUENCE:  CALL  DOUT ( IDUM , IDVEC ) 
where : 

I DUM ....  not  used 
IDVEC  =  .BYTE  -  not  used 

.BYTE  -  module  number,  use  0 
.WORD  -  data  word 

To  set  switches  2  and  6  to  "on",  the  data  word  would  be: 
0000000000100010 

CALLED  BY:  F0R2 . Fortran  program 

PUT. . Fortran  Subroutine 

CALLS: 

DESCRIPTION:  See  program  listing 
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NAME:  GLUE 

LOCATION:  ASSBA . MAC 

DATE:  791112 

AUTHOR:  G.W.M.  Coppus 

ABSTRACT:  Subroutine  takes  two  real  values  and  puts  them  in 
a  four  word  integer  array. 


CALLING  SEQUENCE: 

:  CALL  GLUE(REAL1 ,REAL2, IARRAY) 

where : 


REAL  1 ,  REAL2 

are  two  real  numbers 

IARRAY  is  a  four  word  integer  buffer 


CALLED  BY:  F0R2. 

. . . . For tr an  program 

CALLS: 

« 
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NAME:  GETBC  LOCATION:  SUBF2 . FOR 

DATE:  800218  AUTHOR:  G.W.M.  Coppus 

ABSTRACT:  Subroutine  reads  GC  report  from  HP  1000  computer 
(LARRY)  and  gets  the  value  for  the  bottom  composition  as 
an  integer  number. 

CALLING  SEQUENCE:  CALL  GETBC 

CALLED  BY:  F0R2 . Fortran  program 

CALLS:  SECNDS . System  Subroutine 

MTIN . System  Subroutine 

INDEX . System  Subroutine 

ASCINT . Assembler  Subroutine 

DESCRIPTION: 

First,  the  subroutine  calculates  the  time  available  for 
finding  the  value  of  the  bottom  composition.  This  time 
i  s : 

TREP  =  SAMPLE  TIME  -  GC-DELAY  -  2  =  1 80  -  1 40  -  2  =  38 
seconds 

Then  the  subroutine  starts  looking  for  characters  on  the 
GC  channel,  while  continually  checking  the  elapsed  time, 
unt i 1 : 

1)  The  time  TREP  has  expired.  Then  the  subroutine 
will  exit  whether  it  has  found  a  new  value  for 
the  bottom  composition  or  not.  This  is  done  to 
avoid  a  program  lock  due  to  GC  failure. 

A  carriage  return  character  (octal  15)  is 


2) 


« 


117 


encountered.  Then  the  program  assumes  that  one 
line  of  input  is  complete  and  the  character 
buffer  will  be  searched  for  the  character  string 
&WATER . 

3)  The  character  buffer  is  full,  i.e.  there  are  40 
characters  without  a  carriage  return.  From  the 
GC  report  we  Know  that  the  characterstring 
&WATER  should  occur  within  the  first  40 
characters  of  a  line.  The  character  buffer  will 
be  searched  for  the  character  string. 

If  we  find  the  string  &WATER,  we  Know  that  the  6 
ASCII  characters  preceding  the  string  represent  the 
weight  percentage  of  water  in  the  sample,  e.g. 

94.556  &WATER.  From  this  we  can  deduct  that  the 
methanol  concentration  is  5.444  %wt . 

The  assembler  routine  ASCINT  taKes  care  of  the 
conversion  from  ASCII  to  integer. 

Whether  the  bottom  product  composition  is  found  or 
not,  the  subroutine  will  jump  bacK  in  the  loop  (see 
above)  to  looK  for  new  characters  until  the  time 
TREP  has  expired. 

The  reasons  for  looKing  for  the  water  composition 
and  not  the  methanol  composition  are: 

1)  Sometimes  the  word  " METFIANOL "  is  missing  from 
the  GC  report,  " &WATE R "  seems  to  be  always 
there,  so  this  is  more  reliable. 


ifi 
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2)  Sometimes,  when  the  sample  valve  does  not 

operate  properly,  the  GC  analysis  shows  three 
peaks  instead  of  two.  One  of  the  smaller  peaks 
will  be  identified  as  the  "methanol  peak",  the 
large  peak  is  the  "water  peak".  We  find  that 
100%  -  water  concentration  is  in  this  case  a 
more  accurate  representation  of  the  methanol 
concentration  than  the  number  that  was  found 
from  the  peak  that  was  identified  as  the 
methanol  peak,  so  we  get  increase  accuracy. 

It  should  be  emphasized  that  a  single  erroneous  GC 
reading  will  have  disastrous  results  for  a  run  since 
this  reading  will  cause  a  big  jump  in  the  steam 
flowrate  ,  the  effects  of  which  take  about  one  hour 


to  die  out . 
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NAME:  TIMER  LOCATION:  SUBF2 . FOR 

DATE:  800125  AUTHOR:  G.W.M.  Coppus 

ABSTRACT:  Subroutine  calculates  the  time  for  the  next  sample 
moment  and  suspends  program  execution  until  then. 

CALLING  SEQUENCE:  CALL  TIMER 

CALLED  BY:  F0R2 . Fortran  program 

CALLS:  J ICVT ....  System  Subroutine 

J JCVT ....  System  Subroutine 

JADD . System  Subroutine 

CVTT IM ...  System  Subroutine 
I UNT I L . . . Sys tern  Subroutine 
DESCRIPTION: 

The  beginning  of  the  subroutine  is  straight  forward,  the 
new  sample  time  is  calculated  by  adding  the  sample 
interval  I  SEC  to  the  previous  time.  If,  however,  the 
calculated  time  is  24  hours,  obviously  0  hours 
(midnight)  is  meant,  since  the  internal  clock  of  the 
LSI- 11  switches  to  00:00  after  23:59.  If  for  some  reason 
the  calculated  time  is  later  than  the  present  time,  the 
program  will  continue  immediately,  without  waiting. 
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NAME:  DATAC 


LOCATION:  SUBF2 . FOR 
AUTHOR:  G.W.M.  Coppus 


DATE:  791114 


ABSTRACT:  Subroutine  reads  all  analog  inputs 
CALLING  SEQUENCE:  CALL  DATAC  . 


CALLED  BY:  FOR2 . Fortran  program 


CALLS:  AIHL 


Assembler  Subroutine 


I  ANAL  ....  Assembler  Subroutine 


I PEEKB ...  System  Subroutine 


DESCRIPTION: 

Starting  at  the  first  10  vector  in  the  I/O-table,  the 
subroutine  checks  if  the  vector  is  used  (unused  vectors 
start  with  the  characters  UN). 

If  yes,  the  info  word  in  the  10  vector  (see  section  on 
buffer  usage)  will  be  analyzed.  If  the  vector  represents 
an  analog  input,  the  appropriate  channel  will  be  read. 
Repeat  above  procedure  until  all  vectors  (12)  are  done. 


121 


NAME:  CONTR 


LOCATION:  CONTR. FOR 


DATE:  791114 


AUTHOR:  G.W.M.  Coppus 


ABSTRACT:  Control  algorithm,  calculates  P,  I,  D,  P-Robust 
and  Feed  forward  control  action. 

CALLING  SEQUENCE:  CALL  CONTR ( ERR , SUM , Y ) 
where : 

ERR  =  values  of  error,  YREF-Y 

SUM  =  values  of  error  integral  for  I  action 

Y  =  values  of  controlled  outputs  (top  and  bottom 
compos i t ions ) 

CALLED  BY:  F0R2 . Fortran  program 

CALLS:  ENG . Fortran  function 


Fortran  subroutine 


INGE 


DESCRIPTION: 

This  subroutine  offers  a  variety  of  control  algorithms, 
the  algoritm  is  chosen  by  selecting  ITYPE.  The  meaning 
of  ITYPE  is: 

ITYPE  =  0  no  control  action 


1  proportional  control  2  =  P+feed  forward 


3  Pi-control 


4  =  -+feed  forward 


5  PD-Robust  control 


6  =  -+feed  forward 


7  PID-control 


8  =  -+feed  forward 


9  P-Robust  control 
11  Set  steady  state  values 


10  =  -+feed  forward 
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First  the  program  checks  if  I  TYPE  =  11.  If  yes,  the  last 
output  values  for  feed,  reflux  and  steam  will  replace 
the  old  steady  state  values.  (The  default  values  are: 
Feed  18  g/s,  Reflux  9.5  g/s,  Steam  14  g/s).  Then  I  TYPE 
will  be  set  to  0,  which  means  no  control  (manual).  If 
ITYPE  is  in  the  range  1  -  10,  the  control  subroutine 
will  get  the  feed  flowrate,  top  product  composition  and 
bottom  product  composition  in  engineering  units.  Then 
the  error  vector  will  be  calculated  (YREF-Y).  A 
"computed  GO  TO"  selects  the  next  program  part.  The 
integral  part  and  P-robust  part  have  anti  reset  windup 
which  sets  the  value  of  the  integral  action  to  zero  if 
it  exceeds  a  certain  value.  Feed  forward  action  will  be 
included  for  any  even  value  of  ITYPE. 


■ 
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NAME:  PUT  LOCATION:  SUBF2 . FOR 

DATE:  791114  AUTHOR:  G.W.M.  Coppus 

ABSTRACT:  Subroutine  sends  all  output  signals  to  Current 
Output  Stations. 

CALLING  SEQUENCE:  CALL  PUT 

CALLED  BY:  F0R2 . Fortran  program 

CALLS:  I  ANAL . . . . Assemb 1 er  Subroutine 

DOUT . Assembler  Subroutine 

AOHL . Assembler  Subroutine 

ISLEEP ...  System  Subroutine 
DESCRIPTION: 

This  subroutine  sets  the  Current  Output  Stations  (COS) 
to  the  calculated  values.  There  is  one  analog  output 
connected  from  the  LSI- 11  to  the  COS.  The  proper  COS  is 
selected  by  digital  outputs,  i.e.  closing  switch  number 
1  on  the  digital  output  card  at  the  LSI-11  will  cause 
COS  #1  to  respond  to  the  analog  output  voltage  from  the 
LSI-11. 

The  subroutine  checks  every  10  vector  in  the  I/0-table. 
If  the  vector  is  used  and  an  analog  output,  first  all 
digital  outputs  will  be  set  to  zero  to  blank  out  the 
switching  operation  in  the  LSI,  then  the  analog  output 
will  be  set  to  the  calculated  voltage  and  the 
appropriate  digital  output  will  be  set.  Between  the 
various  switching  operations  in  the  subroutine  there  are 


i  *  I 
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wait  statements  (wait  0.1  secs)  to  enable  the  slow 
swi tches  snd  COS  to  fol low  th©  outputs 
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NAME:  COMP  LOCATION:  SUBF2 . FOR 

DATE:  791112  AUTHOR:  G.W.M.  Coppus 


ABSTRACT:  Completion  subroutine  for  parameter  transmission. 
Sets  a  flag. 

CALLING  SEQUENCE:  parameter  in  CALL  I RCVDF COMP ) 

CALLED  BY:  I RCVDF ...  System  Subroutine 

CALLS: 

DESCRIPTION: 

All  this  subroutine  does  is  setting  FLAGP=.TRUE. 


126 


NAME:  COMD  LOCATION:  SUBF2 . FOR 

DATE:  791112  AUTHOR:  G.W.M.  Coppus 


ABSTRACT:  Completion  subroutine  for  data  transmission.  Sets 
a  flag 

CALLING  SEQUENCE:  parameter  in  CALL  I SDATF COMD ) 

CALLED  BY:  I SDATF ...  System  Subroutine 

CALLS: 

DESCRIPTION: 

All  this  subroutine  does  is  setting  FLAGD=.TRUE. 
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NAME:  ASCINT  LOCATION:  CONV.MAC 

DATE:  791112  AUTHOR:  G.W.M.  Coppus 

ABSTRACT:  Converts  5  ASCII  characters  into  an  integer  number 
CALLING  SEQUENCE:  CALL  ASCINT ( IBOT , IBCOMP ) 
where: 

IBOT  is  a  5  byte  logical  array  containing  5  ASCII  digits 
IBCOMP  is  an  integer  number  0  -  30000 
CALLED  BY:  GETBC . . . . For tran  Subroutine 
CALLS: 

DESCRIPTION: 

This  subroutine  takes  an  array  of  characters 
representing  the  weight  %  of  water  in  a  Gas 
Chromatograph  sample,  e.g.  97555  means  97.555  %wt  water. 
This  number  will  be  converted  to  the  integer  27555.  By 
substracting  this  from  30000  we  get  the  %wt  methanol  in 
the  sample.  An  identical  conversion  (to  integer  97555) 
is  not  possible  because  integer  numbers  only  go  up  to 
32674.  See  description  of  GETBC  for  more  explanation. 


128 


NAME:  I  ANAL  LOCATION:  ASSBA . MAC 

DATE:  791112  AUTHOR:  G.W.M.  Coppus 

ABSTRACT:  Subroutine  analyzes  the  information  bits  and 
channel  number  byte  in  word  4  of  an  10  vector. 

CALLING  SEQUENCE:  CALL  I  ANAL ( IOTAB ( 4 , I ) , 10 , I  AD , ISQ , N ) 
where : 

IOTAB (4, I)  =  info  word  in  10  vector  number  I 
on  return  from  the  subroutine: 

10  =  1  for  an  input 

0  for  an  output 

IAD  =  1  for  analog 
0  for  digital 
ISQ  =  0  for  linear 

1  for  square 

2  for  square  root 

3  for  temperature 

N  =  channel  number  0  -  256 
CALLED  BY:  DATAC . . . . For tran  Subroutine 

ENG . Fortran  Function 

INGE . Fortran  Subroutine 

PUT . Fortran  Subroutine 

CALLS: 

DESCRIPTION: 

Analyzes  the  info  word  in  the  10  vector  bit  by  bit.  See 

of  10  vectors  in  the  section  on  usage  of 


the  description 


buffers. 
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NAME:  AIHL 

LOCATION:  AIHL. MAC 

DATE:  790831 

AUTHOR:  R.L.  Barton 

ABSTRACT:  Subroutine  reads  an  analog  input  on  the  indicated 
channel  and  returns  the  voltage  as  an  integer  number 
representing  millivolts. 

CALLING  SEQUENCE:  CALL  AIHL ( IDUM , IDVEC ) 


where : 

IDUM 

dummy 

IDVEC (  1  )  = 

unused  (set  to  zero) 

IDVEC ( 2 )  = 

channel  #,0-63  on  entry, 

decimal  4  on  return 

IDVEC ( 3  )  = 

number  of  bytes  in  record  (8)  on  return 

I DVEC ( 4 )  = 

address  of  data  buffer  (DATABF)  on  return 

DATABF 

=  .BYTE  count  byte 

.BYTE  Link  Level  Control 

.BYTE  Message  Level  Control 

. BYTE  low  data  byte 

.BYTE  high  data  byte 

.BYTE  not  used 

.BYTE  error  code 

CALLED  BY:  DATAC . . . . For tran  Subroutine 
CALLS: 

DESCRIPTION: 


For  a  detailed  description,  see  the  program  listing.  On 
return  from  the  call,  bytes  4  and  5  in  the  data  buffer 


131 


(DATABF)  contain  the  integer  number  which  represents  the 
voltage  reading.  See  DATAC  to  find  out  how  to  read  from 
the  data  buffer . 


132 


NAME:  ENG  LOCATION:  SUBF2 . FOR 

DATE:  800218  AUTHOR:  G.W.M.  Coppus 

ABSTRACT:  Subroutine  takes  an  10  vector  as  input  and  returns 
the  process  value  in  engineering  units. 

CALLING  SEQUENCE:  VALUE  =  ENG ( IOTAB (1,1) 
where : 

VALUE  =  number  in  engineering  units 
I OT A B (  1  ,  I  )  =  10  vector  from  the  I/0-table. 

CALLED  BY:  CONTR . . . . For tran  Subroutine 
CALLS:  SPLIT  ... .Assembler  Subroutine 

I  ANAL  ....  Assembler  Subroutine 
DESCRIPTION: 

Depending  on  the  type  of  input  (linear  or  square  root  as 
in  flows),  and  the  values  for  SPAN  and  ZERO,  the 
function  will  calculate  the  value  in  engineering  units 
using  the  algorithm: 

IVAL  is  the  integer  representing  voltage  (mVolts) 

RVAL  =  ( IVAL- 1000 ) /4000  for  an  input 
I VAL / 1 0000  for  an  output 
CVAL  =  R VAL**0 . 5 ,  1  or  2  for  SQRT ,  LIN,  SQ 
ENG  =  SPAN*CVAL  +  ZERO 
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NAME:  INGE 

LOCATION:  SUBF2 . FOR 

DATE:  791115 

AUTHOR:  G.W.M.  Coppus 

ABSTRACT:  Calculates  the  integer  value  from  the  number  in 
engineering  units  and  puts  the  number  in  the  10  vector 
Opposi te  of  ENG . 

CALLING  SEQUENCE:  CALL  INGE ( X , IOTAB 1 , I ) 
where : 

X  =  value  in  engineering  units 


IOTAB (1,1)  = 

10  vector  where  integer  is  to  be  put  into 

CALLED  BY:  CONTR , 

. ...Fortran  Subroutine 

CALLS:  SPLIT..., 

.Assembler  Subroutine 

I  ANAL . .  . 

.Assembler  Subroutine 

DESCRIPTION: 

The  algorithm  is: 

CVAL  =  ( X-ZERO ) /SPAN 

RVAL  =  CVAL  for  a  linear  output 

CVAL**0.5  for  a  squared  output 
CVAL**2  for  a  square  root  output 


I VEC ( 3 )  =  1 0000*RVAL 


NAME:  SPLIT 


134 


LOCATION:  ASSBA . MAC 
DATE:  791112  AUTHOR:  G.W.M.  Coppus 


ABSTRACT:  Splits  4  integer  words  into  2  real  values 
CALLING  SEQUENCE:  CALL  SPLIT ( IARRAY , REAL  1 , REAL2 ) 
where : 

IARRAY  =  four  word  integer  array 
REAL1  =  first  real  value 
REAL2  =  second  real  value 

CALLED  BY:  ENG . Fortran  Function 

CALLS: 
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NAME:  AOHL  LOCATION:  AOHL . MAC 

DATE:  791212  AUTHOR:  R.L.  Barton 

ABSTRACT:  Analog  output  driver,  sends  an  analog  output  to 
the  indicated  channel. 

CALLING  SEQUENCE:  CALL  AOHL ( I DUM , I DVEC ) 


where : 

I DVEC  = 

.BYTE  not  used 

.BYTE  module  number  (0) 

.BYTE  converter  on  module  0-3,  we  use  0 

.BYTE  not  used 

.WORD  integer  +10000  -  -10000  (mVolts) 

CALLED  BY:  PUT . Fortran  Subroutine 

CALLS: 

DESCRIPTION: 

See  the  program  listing  for  a  detailed  description.  The 
way  the  control  system  is  set  up,  we  use  module  #0, 
converter  #0  as  the  analog  output  channel  to  drive  the 
Current  Output  Stations. 


ann 
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NAME:  BACK3  LOCATION:  BACK3 . FOR 

DATE:  791112  AUTHOR:  G.W.M.  Coppus 

ABSTRACT:  Mainline  background  program,  the  operator 
interface . 

CALLING  SEQUENCE: 

CALLED  BY: 

CALLS:  ASSIGN  ...  System  Subroutine 
CLOSE ....  System  Subroutine 

GLUE . Assembler  Subroutine 

I RCVDF ...  System  Subroutine 
SPRINT  ...  For tran  Subroutine 
SERVE  ....  For tran  Subroutine 

DESCRIPTION: 

In  the  initialization  section,  all  initial  values  for 
parameters  will  be  read  from  file  INIT.DAT.  Also,  the 
event  file  (EVENT.DAT)  will  be  opened  and  the  time  for 
the  first  scheduled  event  will  be  read. 

The  main  body  of  the  program  consists  of  a  continuous 
loop  around  the  subroutine  SERVE  to  check  for  operator 
commands  and  service  them.  It  also  checks  for  scheduled 
events.  Once  every  sample  time,  the  program  jumps  out  of 
this  loop  to  print  the  just  received  data  on  the  screen. 


. 
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NAME:  SPRINT  LOCATION:  SUBBA . FOR 

DATE.  791110  AUTHOR:  G.W.M.  Coppus 

ABSTRACT:  subroutine  takes  care  of  printing  of  data  on  the 
console 

CALLING  SEQUENCE:  CALL  SPR I  NT ( PCODE ) 
where : 

PCODE  =  0  no  printing 

1  concise  printing 

2  I/O-table 

3  CP-table 

4  PP-table 

5  all  tables 

CALLED  BY:  BACK3 . . . . For tran  program 

SER VE . . . . For tran  Subroutine 

CALLS:  CONC . Fortran  Subroutine 

LI SIO . . . . For tran  Subroutine 
L I SCP . . . . For tran  Subroutine 
L I SPP . . . . For tr an  Subroutine 
DESCRIPTION: 

Depending  on  the  value  of  PCODE,  the  subroutine  calls 
one  or  more  of  the  appropriate  listing  routines. 
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NAME:  CHANPP  LOCATION:  SUBBA . FOR 

DATE:  791112  AUTHOR:  G.W.M.  Coppus 

ABSTRACT:  subroutine  to  enter  changes  in  the  program 
parameters  table. 

CALLING  SEQUENCE:  CALL  CHANPP ( I  PAR , LU ) 
where : 

IPAR  =  2  integers  indicating  the  parameter  to  change 
LU  =  logical  unit  number  to  read  new  values  from 
CALLED  BY:  SERVE  ....  For t ran  Subroutine 
CALLS: 

DESCRIPTION: 

This  subroutine  prompts  the  operator  and  then  reads  new 
parameter  values  into  the  PP-table. 
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NAME:  CHANCP  LOCATION:  SUBBA . FOR 

DATE:  791112  AUTHOR:  G.W.M.  Coppus 

ABSTRACT:  subroutine  to  enter  changes  in  the  controller 
parameters  table. 

CALLING  SEQUENCE:  CALL  CHANCP ( I  PAR , LU ) 
where : 

I  PAR  =  integer  number  of  parameter  to  change 

LU  =  logical  unit  from  where  new  values  are  to  be  read 

CALLED  BY:  SERVE  ....  For t ran  Subroutine 
CALLS: 
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NAME:  CHANIO  LOCATION:  SUBBA . FOR 

DATE:  791112  AUTHOR:  G.W.M.  Coppus 

ABSTRACT:  subroutine  to  enter  changes  in  the  controller 
parameters  table. 

CALLING  SEQUENCE:  CALL  CHAN  10 ( I  PAR , LU ) 
where : 

IPAR  =  two  integers  indicating  the  position  in  the  table 
that  has  to  be  changed. 

LU  =  logical  unit  number  from  where  new  values  are 
to  be  read 

CALLED  BY:  SERVE  ....  For tran  Subroutine 
CALLS:  PUT  IN  ....  For tran  Subroutine 

GLUE . Assembler  Subroutine 


■  ■  •  n 
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NAME:  CONC 

LOCATION:  SUBBA . FOR 

DATE:  791112 

AUTHOR:  G.W.M.  Coppus 

ABSTRACT:  subroutine  lists  a  concise  subset  of  data  on  the 
console 

CALLING  SEQUENCE:  CALL  CONC 

CALLED  BY:  SPR I  NT . . . For t ran  Subroutine 


CALLS:  ENG . 

Fortran  Function 

CVTTIM. . . 

System  Subroutine 

DESCRIPTION: 

This  subroutine  writes  the  most  important  process 
variables  on  the  console  every  sample  time.  At  every 
fifth  sample  a  line  of  mnemonics  will  be  written. 


% 
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NAME:  LISPP  LOCATION:  SUBBA . FOR 

DATE:  791112  AUTHOR:  G.W.M.  Coppus 


ABSTRACT:  subroutine  lists  the  program  parameters  table  on 
the  console 

CALLING  SEQUENCE:  CALL  LISPP 

CALLED  BY:  SPRINT  ...  For tran  Subroutine 

CALLS:  CVTT IM ...  System  Subroutine 


■ 
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NAME:  LISCP  LOCATION:  SUBBA . FOR 

DATE:  791112  AUTHOR:  G.W.M.  Coppus 


ABSTRACT:  subroutine  lists  the  CP-table  on  the  console 

CALLING  SEQUENCE:  CALL  LISCP 

CALLED  BY:  SPRINT  ...  For tran  Subroutine 

CALLS: 
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NAME:  LISIO  LOCATION:  SUBBA . FOR 

DATE:  791112  AUTHOR:  G.W.M.  Coppus 


ABSTRACT:  subroutine  lists  the  I/O-table  on  the  console 
CALLING  SEQUENCE:  CALL  LISIO 
CALLED  BY:  SPR I  NT . . . For tr an  Subroutine 
CALLS:  SPLIT  ....  Assembler  Subroutine 

I  ANAL  ....  Assembler  Subroutine 
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NAME:  PUTIN  LOCATION:  SUBBA . FOR 

DATE:  791112  AUTHOR:  G.W.M.  Coppus 

ABSTRACT:  subroutine  puts  a  new  bit  pattern  in  the  INFO  word 
of  an  10  vector . 

CALLING  SEQUENCE:  CALL  PUTIN(IARG) 
where : 

IARG  =  word  where  bit  pattern  is  to  be  put 
CALLED  BY:  CHANIO . . . For tr an  Subroutine 
CALLS: 

DESCRIPTION: 

This  subroutine  expects  the  operator  to  enter  a  mnemonic 
for  an  10  vector.  From  this  mnemonic  the  subroutine  sets 
up  the  right  bit  pattern.  E.g.  If  the  operator  wants  to 
make  a  vector  an  output  signal  for  a  flow  controller, 
the  mnemonic  would  be  SQRT  AO  (SQuare  Root  Analog  Out) 
and  the  bit  pattern  would  be  00001010.  See  description 
of  buffer  usage. 
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NAME:  SERVE  LOCATION:  SUBBA . FOR 

DATE:  791110  AUTHOR:  G.W.M.  Coppus 

ABSTRACT:  Subroutine  to  service  operator  commands  and 
scheduled  parameter  changes. 

CALLING  SEQUENCE:  CALL  SERVE ( LMPM ) 
where : 

LMPM  =  time  for  next  change  in  minutes  past  midnight 


CALLED 

BY:  BACK3 

. . . .Fortran  program 

CALLS: 

SPRINT. . 

.Fortran  Subroutine 

SCHED. . . 

.Fortran  Subroutine 

CHANPP . . 

.Fortran  Subroutine 

CHANCP . . 

.Fortran  Subroutine 

CHANIO . . 

.Fortran  Subroutine 

I  POKE . . . 

.System  Subroutine 

IPEEK.  .  . 

.System  Subroutine 

ITTINR. . 

.System  Subroutine 

ISDATF. . 

.System  Subroutine 

DESCRIPTION: 

The  operator  can  enter  List  and  Change  commands  on  the 
console.  The  subroutine  analyzes  the  command  and  invokes 
the  appropriate  subroutine.  Also,  the  subroutine  checks 
the  time  of  scheduled  events,  and  reads  the  event  from 
disk  file  at  the  right  time  for  implementation. 
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NAME:  SCHED  LOCATION:  SCHED . FOR 

DATE:  791218  AUTHOR:  G.W.M.  Coppus 

ABSTRACT:  subroutine  reads  and  implements  events  from  an 
event  file. 

CALLING  SEQUENCE:  CALL  SCHED ( LMPM , SEND ) 
where : 

LMPM  =  time  for  next  event  in  minutes  past  midnight 
SEND  =  logical  flag  to  indicate  a  change 
CALLED  BY:  SERVE  ....  For tran  Subroutine 
CALLS:  CVTT IM ...  System  Subroutine 
CHANIO ...  For tran  Subroutine 
CHANCP ...  For tran  Subroutine 
CHANPP ...  For tran  Subroutine 
DESCRIPTION: 

The  subroutine  checks  the  next  sample  time  and  compares 
it  with  the  time  for  the  next  change.  If  the  time  for 
the  next  change  is  later  than  the  next  sample  time,  or 
not  within  ten  minutes  before  the  next  sample  time, 
nothing  will  happen. 

Otherwise,  the  new  parameters  will  be  read  from  the 
event  file  and  subroutine  SERVE  will  be  flagged  to  send 
them  to  the  foreground  program. 


7.  APPENDIX  B:  Plots  of  simulation  and  experimental 

responses 
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Plot  1.  Open  loop,  +5%  steam  flow  rate 

SIMULATION  RUN  (flows  in  grams/sec,  compositions 
weight  fraction  methanol) 


in 


CLOSED  LOOP  SlflULRTION  CLOSED  LOOP  SIMULATION 


150 


© 

o 


l/J 

© 


o 


l/> 

M 


o 

o 


U5  M 
^  h~ 


o 

u> 


W 


o 


ncru  uw3is 


2 


a 

cH 


a 

_J 

o 


-dU03  1-0X108 


2 

o 


cr 

_J 

ro 

n 


a. 

o 

o 


a 

UJ 

to 

o 


Plot  2. 


Open  loop,  "5% 
SIMULATION  RUN 
weight  fraction 


steam  flow  rate 

(flows  in  grams/sec,  compositions 
methanol ) 
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Plot  3 


P-control ,  +5%  feed  flow  rate 

SIMULATION  RUN  (flows  in  grams/sec,  compositions  in 
weight  fraction  methanol) 
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Plot  4.  P-control ,  -5%  feed  flow  rate 

'  SIMULATION  RUN  (flows  in  grams/sec,  compositions  in 
weight  fraction  methanol) 


CLOSED  LOOP  SIMULATION 
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Plot  5.  P-robust  feedback  control,  -20%  feed  flow  rate 

SIMULATION  RUN  (flows  in  grams/sec,  compositions  in 
weight  fraction  methanol) 
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Plot  6.  PI  control,  -20%  feed  flow  rate 

SIMULATION  RUN  (flows  in  grams/sec,  compositions  in 
weight  fraction  methanol) 
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Plot  7.  P-robust  feedback  contro 
compos i t i on 

SIMULATION  RUN  (flows  in 
weight  fraction  methanol 
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Plot  8.  PI  control,  +1%  top  product  composition 

SIMULATION  RUN  (flows  in  grams/sec,  compositions  fin 
weight  fraction  methanol) 
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Plot  9.  P-robust  feedback  control,  +1%  bottom  product 
compos i t i on 

c  T  ml  J I  ATION  RUN  (flows  in  grams/sec,  compositions  in 
weight  fraction  methanol) 
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Plot  10.  PI  control,  +1%  bottom  product  composition 

SIMULATION  RUN  (flows  in  grams/sec,  compositions  in 
weight  fraction  methanol) 
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Plot  11.  P-robust  feedback  feedforward  control,  -20%  feed 
flow  rate 

SIMULATION  RUN  (flows  in  grams/sec,  compositions  in 
weight  fraction  methanol) 
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Plot  12.  PI -plus-feedforward  control,  -20%  feed  flow  rate 

SIMULATION  RUN  (flows  in  grams/sec,  compositions  in 
weight  fraction  methanol) 
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Plot  13.  PD-robust  feedback  feedforward  control,  -20%  feed 
flow  rate 

SIMULATION  RUN  (flows  in  grams/sec,  compositions  in 
weight  fraction  methanol) 
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Plot  14.  PID-plus-feedforward  control,  -20%  feed  flow  rate 

SIMULATION  RUN  (flows  in  grams/sec,  compositions  in 
weight  fraction  methanol) 
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Plot  15.  PD-robust  feedback  feedforward  control,  +1%  top 
product  composition 

SIMULATION  RUN  (flows  in  grams/sec,  compositions  in 
weight  fraction  methanol) 
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Plot  16.  P I D - p 1  us- feedforward  control,  +1%  top  product 
compos i t i on 

SIMULATION  RUN  (flows  in  grams/sec,  compositions  in 
weight  fraction  methanol) 
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Plot  17.  PD-robust  feedback  feedforward  control,  +1%  bottom 
product  composition 

SIMULATION  RUN  (flows  in  grams/sec,  compositions  in 
weight  fraction  methanol) 


CLOSED  LOOP  SIMULATION  '  CLOSED  LOOP  SIMULATION 


166 


o 


o 


a: 

-J 

13 

n 

H 

CO 

CL 

o 

o 


a 

0 

o 

_j 

u 


o 

►  -« 
>- 
cr . 

ID 

n 

M 

CO 

a. 

a 

o 


n 

UJ 

8 

d 


raid  Tutmnr.iu 


'tJlJOD  dOi 


Plot  18.  PID-plus-feedforwand  control,  +1%  bottom  product 
compos i t i on 

SIMULATION  RUN  (flows  in  grams/sec,  compositions 
weight  fraction  methanol) 
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Plot  19.  Open  loop,  +10%  reflux  flow  rate 

EXPERIMENTAL  RUN  (flows  in  grams/sec,  compositions 
in  weight  %  methanol ) 
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Plot  20.  Open  loop,  -10%  reflux  flow  rate 

EXPERIMENTAL  RUN  (flows  in  grams/sec,  compositions 
in  weight  %  methanol) 
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Plot  21.  P  control  (gains  3.571  top  and  -0.576  bottom),  +20% 
feed  flow  rate 

EXPERIMENTAL  RUN  (flows  in  grams/sec,  compositions 
in  weight  %  methanol) 
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Plot  22.  P  control  (gains  2‘.  5  top  apd  -0.4  bottom),  +10% 
feed  flow  rate 

EXPERIMENTAL  RUN  (flows  in  grams/sec,  compositions 
in  weight  %  methanol ) 
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23.  P*  control  (gains  2.5  top  and  -0.4  bottom),  -10% 
feed  flow  rate 

EXPERIMENTAL  RUN  (flows  in  grams/sec,  compositions 
in  weight  %  methanol) 
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Plot  24. 


P  control  (gains  1.786  top 
feed  flow  rate 
EXPERIMENTAL  RUN  (flows  in 
in  weight  %  methanol) 


and  -0.288  bottom),  +10% 
grams /sec,  compositions 
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25.  P-robust  feedback  feedforward  control,  -25%  feed 
f low  rate  . 

EXPERIMENTAL  RUN  (flows  in  grams/sec,  compositions 
in  weight  %  methanol) 
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Plot  26.  P-robust  feedback  feedforward  control,  +25%  feed 
f low  rate  ,  . 

EXPERIMENTAL  RUN  (flows  in  grams/sec,  compositions 
in  weight  %  methanol) 
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Plot  27.  p-robust  feedback  feedforward  control,  +1%  top 
product  composition 

EXPERIMENTAL  RUN  (flows  in  grams/sec,  compositions 
in  weight  %  methanol) 
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Plot  28. 


P-robust  feedback  feedforward  control,  -1%  top 
product  composition 

EXPERIMENTAL  RUN  (flows  in  grams/sec,  compositions 
in  weight  %  methanol) 
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Plot  29. 
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P-robust  feedback  feedforward  control,  +1%  bottom 
product  composition 

EXPERIMENTAL  RUN  (flows  in  grams/sec,  compositions 
in  weight  %  methanol ) 
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Plot  30.  R-robust  feedback  feedforward  control,  -1%  bottom 
product  composition 

EXPERIMENTAL  RUN  (flows  in  grams/sec,  compositions 
in  weight  %  methanol) 
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Plot  31.  PI  control,  -25%  feed  flow  rate 

EXPERIMENTAL  RUN  (flows  in  grams/sec,  compositions 
in  weight  %  methanol ) 
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Plot  32.  PID  control,  -25%  feed  flow  rate 

EXPERIMENTAL  RUN  (flows  in  grams/sec,  compositions 

in  weight  %  methanol) 
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Plot  33.  PID-plus-feedforward  control,  -25%  feed  flow  rate 
EXPERIMENTAL  RUN  (flows  in  grams/sec,  compositions 
in  weight  %  methanol) 
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Plot  34.  P ID-plus- feedforward  control,  +1%  top  product 
compos i t i on 

EXPERIMENTAL  RUN  (flows  in  grams/sec,  compositions 
in  weight  %  methanol) 
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35.  P I D-p 1  us- feedforward  control,  +1%  bottom  product 
compos i t i on 

EXPERIMENTAL  RUN  (flows  in  grams/sec,  compositions 
in  weight  %  methanol) 


Plot  36.  PD-robust  feedback  feedforward  control,  +25%  feed 
flow  rate 

EXPERIMENTAL  RUN  (flows  in  grams/sec, 
in  weight  %  methanol) 
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37.  PD-robust  feedback  feedforward  control,  -25%  feed 
flow  rate 

EXPERIMENTAL  RUN  (flows  in  grams/sec,  compositions 
in  weight  %  methanol) 
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Plot  38 
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PD-robust  feedback  feedforward  control,  +1%  top 
product  composition 

EXPERIMENTAL  RUN  (flows  in  grams/sec,  compositions 
in  weight  %  methanol) 


187 


cUDD  HJU 


r~ 
88  iZ 


- r~ 

oa  a z 


£0  91 
*<nj  <naj 


ea  zi 


8  S 


< 


- 1 - - - r 

EO  96  00  16  08  C6 

<303  dUi 


Plot  39? 


PD-robust  feedback  feedforward  control,  -1%  top 
product  composition 

EXPERIMENTAL  RUN  (flows  in  grams/sec,  compositions 
in  weight  %  methanol) 
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Plot  40. 


PD-robust  feedback  feedforward  control,  +1%  bottom 
product  composition 

EXPERIMENTAL  RUN  (flows  in  grams/sec,  compositions 
in  weight  %  methanol) 
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41.  PQ-robust  feedback  feedforward  control,  -1%  bottom 
product  composition 

EXPERIMENTAL  RUN  (flows  in  grams/sec,  compositions 
in  weight  %  methanol) 
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